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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows: 


Multiply 


By 


To  Obtain 


cubic  feet 


0.02831685 


cubic  meters 


feet 


0.3048 


meters 


miles  (U.S.  statute) 


1.609347 


kilometers 


1  Introduction 


Background 

The  U.S.  Army  Corps  of  Engineers  (CE)  controls,  maintains,  and  conserves 
water  resources  to  provide  for  flood  control,  navigation,  irrigation,  power  gen¬ 
eration,  recreation,  water  quality,  water  supply,  fish  and  wildlife  protection  and 
enhancement,  and  other  authorized  project  purposes  by  regulating  releases 
from  dams.  The  degree  to  which  preimpoundment  flows  are  regulated  depends 
upon  project  purposes.  Alteration  of  preimpoundment  flows  ranges  from  rela¬ 
tively  long-term  seasonal  changes  required  to  support  navigation  or  to  provide 
for  downstream  flood  protection  to  extremely  short-term  changes  required  for 
peaking  hydropower  generation. 

Increased  economic  development  in  river  basins,  plus  changing  demogra¬ 
phic,  social,  and  land  use  patterns,  is  placing  increased  demands  on  the  natural 
resources  in  river  systems.  Operation  of  dams  impacts  several  natural  resource 
categories  including  inpool  reservoir  fish  resources,  wetlands,  wildlife,  and 
downstream  fish  resources.  Efforts  by  the  CE  to  continue  fostering  economic 
development  in  river  basins,  while  simultaneously  protecting  environmental 
quality,  require  the  availability  of  predictive  tools  Aat  can  be  used  to  balance 
the  developmental  and  environmental  needs  of  the  region. 

This  report  documents  the  approaches  and  underlying  rationale  used  to 
comparatively  assess  the  effects  of  physical  habitat  changes  resulting  from 
alternative  reservoir  operations  on  the  native  warmwater  riverine  fish  commu¬ 
nity  in  the  tailwaters  of  the  main  stem  dams  on  the  Missouri  River.  The  main 
stem  Missouri  River  reservoir  system  extends  from  Montana  to  Nebraska  and 
includes  the  following  reservoir  projects:  Fort  Peck  Dam  (River  Mile  (RM) 
1771.5),  Garrison  Dam  (RM  1389.9),  Oahe  Dam  (RM  1072.3),  Big  Bend  Dam 
(RM  987.4),  Fort  Randall  Dam  (880.0),  and  Gavins  Point  Dam  (RM  811.1). 
However,  results  of  the  analysis  are  presented  for  the  tailwater  of  Gavins  Point 
Dam  only  because  it  was  the  most  thoroughly  studied  and  documented  tail- 
water.  Analysis  of  the  other  reservoirs  on  the  main  stem  Missouri  River  hav¬ 
ing  a  downstream  tailwater  (instead  of  releasing  into  the  headwaters  of  a 
downstream  reservoir)  can  be  foimd  in  Nestler,  Schneider,  and  Hall  (1993). 
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Problem 


Like  many  CE  Districts  and  Divisions,  the  Missouri  River  Division  (MRD) 
lacked  systematic  assessment  tools  necessary  to  analyze  the  effects  of  reservoir 
operation  on  physical  habitat  for  native  riverine  fishes  occurring  in  river 
reaches  affected  by  reservoir  operation.  The  lack  of  assessment  tools  pre¬ 
vented  MRD  from  developing  and  documenting  operational  plans  that  optimize 
economic  and  natural  resource  categories  in  the  main  stem  Missouri  River. 
Assessment  tools  are  required  to  predict  and  manage  the  environmental  effects 
of  reservoir  operation  to  facilitate  optimal  use  of  this  valuable  resource.  For 
example,  a  major  water  diversion  may  potentially  impact  fish  in  the  system  by 
altering  flow  regimes  and  water  quality.  Conversely,  implementation  of  mea¬ 
sures  beyond  those  required  to  maintain  and  protect  environmental  quality  of 
the  system  may  preclude  drawing  benefits  from  other  uses  of  water  in  the  main 
stem  system. 


Objective 

The  objective  of  this  study  was  to  develop  a  method  that  could  be  used  to 
develop  simple  relationships  between  flow  and  habitat  values  allowing  CE 
engineers  and  scientists  to  predict  the  tailwater  habitat  value  downstream  of 
reservoirs  for  different  operational  alternatives  or  channel  configurations.  The 
method  was  specifically  designed  to  describe  the  effects  of  reservoir  operation 
on  tailwater  fish  habitat  in  a  context  suitable  for  reservoir-specific  or  basin¬ 
wide  trade-off  analysis.  The  method  developed,  the  Riverine  Community 
Habitat  Assessment  and  Restoration  Concept  (RCHARC),  is  applied  to  the 
Gavins  Point  Dam  tailwater  as  a  case  history  that  can  be  emulated  for  analysis 
of  other  reservoirs  that  release  steady  or  gradually  varied  flows.  Modifications 
to  the  RCHARC  to  allow  assessment  of  peaking  releases  associated  with 
hydropower  operation  can  be  found  in  Nestler,  Schneider,  and  Hall  (1993). 


Model  Strategy 

Analyzing  environmental  effects  requires  long-term  time  series  estimates  of 
the  releases  from  the  dam  associated  with  each  operational  alternative.  Usually 
the  discharge  estimates  are  obtained  from  a  reservoir  operation  or  basin  hydrol¬ 
ogy  model.  Each  alternative  operation  is  represented  as  a  set  of  rules  that 
determines  the  operation  of  the  dam  as  a  function  of  reservoir  stages  or  down¬ 
stream  flow  requirements.  The  set  of  rules  is  applied  to  a  hydrologic  time 
series  of  inflows  into  the  reservoir  or  reservoir  system  to  determine  the  pattern 
of  releases  associated  with  each  alternative.  The  summary  environment^ 
effects  of  each  alternative  are  determined  by  integrating  the  environmental 
effects  of  each  release  at  a  particular  time-step  (usually  1  month)  over  the 
hydroperiod. 
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For  this  application  of  the  RCHARC,  MRD  used  the  Long-Range  Study 
(LRS)  model  to  predict  monthly  reservoir  stages  and  monthly  reservoir  dis¬ 
charges  associated  with  each  operational  alternative  (U.S.  Army  Engineer 
Division,  Missouri  River,  1992).  The  operating  rules  for  each  alternative  were 
each  run  separately  by  the  LRS  model  using  the  93-year  hydrological  record 
available  for  the  Missouri  River.  The  LRS  model  generated  a  unique  set  of 
stages  and  discharges  characteristic  of  each  alternative  that  were  the  basis  of 
the  habitat  analysis. 

Use  of  a  hydrologic  model  to  generate  stage  and  discharge  information  for 
each  reservoir  ordinarily  limits  the  data  available  to  perform  the  habitat  analy¬ 
sis.  The  only  data  available  to  evaluate  the  alternatives  are  the  unique  monthly 
stages  and  discharges  available  for  each  alternative  and  data  that  are  associated 
with  the  93-year  hydrological  record.  Time  series  of  measured  meteorological 
and  water  quality  data  are  usually  not  as  complete  as  the  time  series  of  stages 
and  discharges  predicted  by  the  hydrology  models.  Consequently,  the  app¬ 
roach  used  to  model  tailwater  habitat  was  restricted  to  those  variables  that 
were  either  predicted  by  the  LRS  model  or  whose  trends  could  be  derived 
from  historical  data  linked  to  the  93-year  hydrologic  record. 

For  this  application  of  the  RCHARC,  the  software  estimating  the  habitat 
impact  of  each  alternative  was  coded  as  a  separate  module  that  was  linked  to 
the  LRS  model.  The  LRS  model  was  run  once  for  each  alternative  over  the 
93-year  hydroperiod,  and  the  output  was  passed  to  the  module  deteimining  the 
monthly  reservoir-specific  tailwater  habitat  values  for  each  flow  alternative. 
Other  resource  categories  can  be  assessed  in  a  conceptually  parallel  method; 
i.e.,  each  resource  category  is  assessed  by  a  separate  module  that  is  fed  infor¬ 
mation  from  a  common  reservoir  operation  or  basin  hydrology  model. 

It  is  often  necessary  or  useful  that  the  results  of  a  habitat  analysis  done  as 
part  of  a  trade-off  assessment  be  expressed  as  a  “value  function;”  i.e.,  each 
different  average  discharge  has  an  associated  habitat  value  ranging  from  0.0  to 
1.0  with  0.0  representing  minimum  value  and  1.0  representing  maximum 
value.  Expressing  the  relationships  between  flow  and  different  resource  cate¬ 
gories  as  value  functions  allows  direct  comparison  of  seemingly  disparate 
resource  categories  such  as  recreation,  warmwater  fish  habitat,  or  flood  control 
benefits  that  are  assessed  using  different  units  of  measurements. 

The  habitat  assessment  of  the  Gavins  Point  tailwater  using  the  RCHARC 
has  the  following  steps: 

a.  The  LRS  model  predicts  monthly  flow  estimates  for  each  alternative 
over  the  93-year  period  of  record. 

b.  Separate  value  functions  are  formulated  for  each  month. 

c.  Each  monthly  flow  predicted  by  the  LRS  model  is  replaced  by  an  asso¬ 
ciated  habitat  value. 
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d.  The  effect  of  each  alternative  is  obtained  by  summing  (or  integrating  in 
some  other  way,  e.g.,  using  annual  habitat  minimums  in  cases  where 
the  armual  flow  pattern  differs  substantially  from  the  comparison  stan¬ 
dard)  the  resulting  habitat  values  over  the  93-year  period  of  record. 

e.  The  alternatives  are  evaluated  against  each  other  or  against  a  standard 
to  rank  the  alternatives  and  thus  provide  the  rationale  for  selecting 
optimal  operational  plans  as  part  of  a  trade-off  analysis. 

Assessing  the  effects  of  reservoir  operation  can  be  complicated  by  the 
choice  of  the  “no-action”  alternative.  For  this  application  of  the  RCHARC, 
MRD  determined  that  the  no-action  alternative  is  defined  by  the  present  system 
of  rules  currently  used  to  operate  the  main  stem  Missouri  River  dams.  Ordi¬ 
narily,  the  no-action  alternative  serves  as  an  implicit  standard  against  which 
other  alternatives  can  be  evaluated.  In  this  case,  the  no-action  alternative  is 
itself  an  operational  alternative  (and  not  an  unregulated  condition)  whose  habi¬ 
tat  value  has  not  been  quantified.  The  no-action  alternative  was  treated  as  one 
of  the  operational  alternatives,  and  its  habitat  impact  was  evaluated  by  compar¬ 
ing  its  depth/velocity  distributions  to  the  “comparison  standard.”  The  relative 
impacts  of  alternative  operations  were  then  evaluated  relative  to  the  score  of 
the  no-action  alternative.  Thus,  the  habitat  impact  of  each  alternative,  includ¬ 
ing  the  no-action  alternative,  was  first  quantified  against  the  “standard  of  com¬ 
parison.”  The  habitat  value  of  each  alternative  flow,  including  the  no-action 
alternative,  was  then  ranked  using  the  standard  of  comparison  as  the  impact 
gauge.  The  habitat  impacts  of  the  project  alternatives  were  evaluated  against 
the  habitat  value  of  the  no-action  alternative  to  complete  the  analysis. 


Approach 

The  “value  function”  assessment  framework  developed  by  the  MRD 
requires  that  a  systematic  procedure  be  developed  to  determine  the  habitat 
value  of  each  monthly  flow  for  the  tailwater.  Various  methods  are  available  to 
relate  flow  to  fish  habitat  value  (Reiser,  Wesche,  and  Estes  1989).  The  most 
commonly  applied  methodology,  the  Physical  Habitat  Simulation  System 
(PHABSIM)  of  the  Instream  Flow  Incremental  Methodology  (IFIM),  developed 
by  the  U.S.  Fish  and  Wildlife  Service  (USFWS)  (Milhous,  Updike,  and 
Schneider  1989),  has  received  considerable  criticism  for  applications  on  large, 
warmwater  river  systems  (Bain  and  Boltz  1989;  Nestler  1993).  Even  if  the 
IFIM  were  considered  defensible  for  this  application,  the  short  time  frame 
available  for  this  study  prevents  its  application  for  the  numerous  warmwater 
fishes  occurring  in  the  Missouri  River  system. 

Suitability  curves  relate  the  habitat  value  of  different  depths  and  velocities 
to  species  life  stages  or  guilds  (Bovee  1986).  Major  alterations  in  channel 
morphology  and  flow  regime  make  development  of  defensible  suitability 
curves  problematic.  As  Tyus  (1992)  points  out,  it  is  difficult  to  develop 
suitability  curves  for  impacted  species  in  highly  modified  systems  because  the 
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system  is  deficient  in  one  or  more  critical  habitat  components  or  else  the 
impacted  species  would  not  be  rare.  It  is  difficult  or  impossible  to  develop 
suitability  curves  for  habitat  features  that  are  rare  or  missing.  Therefore,  it  is 
unlikely  that  habitat  suitability  curves  developed  from  a  highly  impacted  sys¬ 
tem  like  the  Missouri  River  will  capture  all  necessary  habitat  features  required 
by  the  target  species. 

Additionally,  developing  a  single  value  that  relates  the  value  of  a  particular 
flow  alternative  to  the  entire  warmwater  fish  community  would  be  nearly 
impossible  using  suitability  curves  as  the  basis  of  the  analysis.  The  approxi¬ 
mately  60  species  and  four  life  stages  per  species  require  the  flow  require¬ 
ments  of  nearly  250  species’  life  stages  be  determined,  summarized,  and 
integrated.  Creation  of  defensible  suitability  curves  for  even  one  species  can 
be  difficult  and  time-consuming  (Bovee  1986).  Formulating  habitat  suitability 
information  on  even  a  few  species  in  large  rivers  can  be  a  particularly  expen¬ 
sive  and  protracted  activity  because  of  the  difficulty  in  sampling  at  the  large 
temporal  and  spatial  scales  characteristic  of  large  river  systems. 

While  use  of  suitability  information  to  define  the  habitat  requirements  of 
large  warmwater  fishes  does  not  seem  to  be  technically  desirable  or  logistically 
feasible,  sufficient  literature  is  available  that  suggests  that  fishes  respond  to 
depth  and  velocity  patterns  including  fishes  in  large  warmwater  rivers  (Bain, 
Reed,  and  Scheidegger  1991).  Bain,  Reed,  and  Scheidegger  (1991)  point  out 
that  the  distribution  of  warmwater  fishes  in  the  Cahaba  River  system  of 
Alabama,  a  large,  southeastern,  warmwater  river  system  supporting  fish  species 
that  also  occur  in  the  Missouri  River,  can  be  largely  explained  using  broad 
depth  and  velocity  categories.  The  work  of  Bain,  Reed,  and  Scheidegger  sug¬ 
gests  that  broad  depth  and  velocity  information  should  be  included  in  a  habitat 
analysis.  However,  because  of  the  limitations  associated  with  the  use  of  habi¬ 
tat  suitability  information,  depth/velocity  information  should  be  incorporated  in 
a  marmer  that  does  not  use  species  life  stage  or  guild  habitat  suitability  curves. 

A  large  warmwater  river  system  such  as  the  Cahaba  River,  Alabama,  can 
have  nearly  100  species,  each  of  which  can  have  several  life  stages.  Attempt¬ 
ing  to  create  the  necessary  suitability  curves,  determine  the  flow-habitat  rela¬ 
tionships  based  on  each  curve,  integrating  the  population  level  information  to 
make  community-level  statements,  and  then  distilling  this  information  to  facili¬ 
tate  wise  decision  making  is  a  daunting  task.  A  more  feasible  approach  is  to 
find  a  reach  of  river  that  contains  a  healthy  community.  The  physical  habitat 
patterns  over  an  armual  cycle  of  this  river  reach  can  then  be  used  as  a  standard 
to  which  the  habitat  conditions  in  an  impacted  system  can  be  compared. 
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2  Riverine  Community 
Habitat  Assessment  and 
Restoration  Concept 


Based  on  considerations  discussed  in  the  previous  section,  RCHARC,  a  new 
approach  for  relating  the  effects  of  flow  alterations  and  alternative  channel 
designs  on  aquatic  biota,  was  developed.  This  system  combines  conceptual 
elements  of  the  Index  of  Biotic  Integrity  GBI)  (Karr  et  al.  1986)  and  the 
PHABSIM  system.  Like  the  IBI,  the  RCHARC  requires  use  of  a  river  system 
as  toe  basis  of  comparison,  that  is,  a  “comparison  standard”  for  toe  analysis 
against  which  toe  various  project  alternatives  can  be  evaluated.  The  compari¬ 
son  standard  river  system  (CSRS)  is  considered  to  represent  toe  ideal  habitat 
conditions,  both  in  terms  of  channel  configuration  and  seasonally  varying  flow 
characteristics,  for  toe  aquatic  community  in  toe  project  river  system.  The 
CSRS  can  be  selected  based  on  professional  consensus,  physical  similarity  to 
toe  project  system,  or  similarity  of  toe  aquatic  community  in  toe  standard 
system  to  what  is  desired  in  toe  project  system.  The  CSRS  can  be  a  nearby 
river  system,  reaches  of  the  river  further  upstream  or  downstream  of  toe  proj¬ 
ect  and  not  impacted  by  toe  project,  or  toe  project  river  reach  but  evaluated  in 
a  “witoout-project”  condition. 

For  this  study,  toe  CSRS  was  considered  to  be  toe  preproject  Missouri 
River  channel  under  preproject  flow  conditions.  This  choice  was  based  on 
several  considerations.  First,  channel  modification  resulting  from  flow  regula¬ 
tion  and  dismption  of  downstream  sediment  transport  has  been  documented  on 
numerous  occasions  (e.g.,  U.S.  Army  Engineer  District,  Omaha,  1989).  There 
seems  to  be  concomitant  decline  in  those  species  that  comprise  toe  native 
warm  water  fish  community  of  this  river  (Gardner  and  Stewart  1987;  Pfleiger 
and  Grace  1987;  Schmulbach  1974;  and  Hesse  et  al.  1989).  These  same 
authors  have  speculated  that  toe  decline  in  these  species  is  related  to  the 
impacts  of  regulation  on  a  variety  of  factors,  including  alteration  in  down¬ 
stream  physical  habitat. 

This  study  assumes  that  a  return  to  preproject  habitat  conditions  for  both 
flow  and  channel  characteristics  should  provide  optimum  physical  habitat  for 
the  warmwater  fish  community  of  toe  Missouri  River,  all  other  factors  being 
equal.  Many  ecological  processes  are  influenced  by  depth  and  velocity 
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conditions  in  the  channel,  such  as  organic  matter  transport,  substrate  composi¬ 
tion,  and  bed  form,  so  that  a  return  to  preproject  flow  conditions  should  also 
tend  to  return  these  flow-associated  factors  closer  to  their  preproject  conditions. 
In  the  context  of  an  RCHARC  assessment,  a  project  alternative  that  provides 
depth,  velocity,  and  top  width  conditions  closer  to  the  CSRS,  in  this  case  the 
preproject  Missouri  River,  will  rank  higher  than  a  project  alternative  that  pro¬ 
vides  less  similar  conditions. 

Although  the  RCHARC  does  not  directly  assess  impacts  on  species  using 
life  stage-specific  suitability  curves,  this  information  is  implicit  in  the  method¬ 
ology.  For  a  given  flow,  a  distribution  of  depth  and  velocity  is  associated  with 
that  flow  (Figure  1,  left  ordinate).  These  distributions  represent  the  habitat 
template  upon  which  the  community  is  structured.  As  Bain,  Reed,  and 
Scheidegger  (1991)  have  pointed  out,  some  species  will  ordinate  along  shallow 
portions  of  the  depth  distribution  and  others  will  ordinate  out  on  deeper  por¬ 
tions  of  the  depth  distribution.  Thus,  the  composition  of  the  fish  community 
will  be  determined  by  the  distribution  of  depth/velocity,  all  other  factors  being 
equal.  Changes  in  the  frequency  distribution  of  depth  and  velocity  wiU  result 
in  associated  changes  in  the  warmwater  fish  community.  A  shift  in  the  fre¬ 
quency  distribution  of  depth  that  reduces  the  amount  of  shallow  water  will 
favor  species  that  require  deeper  water  (Figure  1). 

Competition  and  predation  are  included  indirectly  in  RCHARC  since  the 
physical  habitat  of  the  standard  system  also  provides  the  template  upon  which 
interspecies  interactions  such  as  competition  and  predation  are  structured. 

Thus,  duplication  of  the  physical  structure  of  the  standard  system  in  the  project 
system  wiU  also  provide  a  template  conducive  to  the  duplication  of  species 
interactions. 

The  RCHARC  summarizes  physical  habitat  in  rivers  fundamentaUy  differ¬ 
ently  from  the  methods  employed  in  the  PHABSIM  system.  NoimaUy, 
instream  flow  studies  using  a  habitat-based  approach  (the  IFIM  representing 
the  most  common  approach)  use  point-by-point  measures  or  predictions  of 
depth  and  velocity  in  ceUs  of  cross  sections  under  different  discharges  and 
relate  each  point  to  its  suitability  to  a  particular  target  life  stage.  The  suitabU- 
ity  curve  in  turn  parallels  the  idea  of  simulating  and  evaluating  habitat  on  a 
point-by-point  basis.  That  is,  each  point  is  evaluated  relative  to  a  scale  of 
0  to  1  without  consideration  of  other  points  (Figure  2).  However,  biologists 
have  often  observed  that  the  habitat  requirements  for  many  aquatic  species  are 
based  on  more  than  simple  point-by-point  evaluations.  For  example,  a  preda¬ 
tory  fish  may  select  deeper,  slower  water  as  a  velocity  shelter,  but  the  slower 
water  wiU  often  be  near  faster  water  that  wiU  transport  drift  or  other  food 
entrained  in  the  water  (Hughes  and  Dill  1990).  A  more  complete  analysis  may 
indicate  that  much  of  the  food  eaten  by  the  predatory  fish  may  be  produced  in 
a  specific  part  of  the  river  such  as  a  riffle. 

Depiction  of  the  habitat  requirements  of  stream  fishes  in  a  methodology 
that  focuses  on  point-by-point  analyses  is  unrealistic  because  there  is  no 
explicit  mechanism  to  include  habitat  diversity  in  the  assessment. 
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or  velocity  increment.  The  right  ordinate  represents  the  habitat  value  from  0.0  to  1.0  for  each  depth  or  velocity.  The  relative 
value  of  the  habitat  for  each  species  can  be  determined  by  how  much  of  the  frequency  distribution  falls  within  its  suitability 
curve,  or  restated,  the  composition  of  this  hypothetical  group  of  species  is  at  least  partially  determined  by  the  distribution  of 
depth  and  velocity 
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Figure  2.  Comparison  of  a  single  population  or  guild  analysis  with  the  RCHARC  approach. 

Both  approaches  begin  with  channel  cross-section  information  (cross  section  at 
three  discharges,  Q1 ,  Q2,  and  Q3)  with  cell-specific  depth  and  velocity  information 
available  in  each  cell  at  the  different  discharges.  In  the  usual  habitat-based 
approach,  the  cell-by-cell  depths  and  velocities  at  each  discharge  are  separately 
evaluated  against  suitability  curve  information.  The  habitat  value  of  a  particular 
cross  section  at  a  particular  discharge  is  obtained  by  summing  the  habitat  values  of 
its  component  cells  to  generate  the  plot  in  the  upper  right.  The  plot  relates  the 
habitat  value  of  the  transect  at  each  discharge  to  the  habitat  requirements  of  the 
biological  target.  However,  expanding  this  approach  to  a  community  or  ecosystem 
level  can  be  difficult  and  cumbersome.  In  the  RCHARC  approach,  the  cell-by-cell 
depth  and  velocity  information  for  each  discharge  is  summarized  with  a  frequency 
distribution.  The  three  resulting  frequency  distributions  are  then  compared  to  the 
frequency  distribution  of  the  habitat  standard  to  determine  how  the  alternatives  rank 
to  the  standard  (plot  in  lower  right) 

Point-by-point  analysis  may  show  that  a  stream  with  relatively  uniform 
hydraulic  conditions  (i.e.,  a  ditch)  provides  optimum  habitat  for  a  particular 
life-stage  because  all  of  the  depth  and  velocity  conditions  fall  within  the  maxi¬ 
mum  habitat  value  of  the  suitability  curve.  Conversely,  a  stream  containing 
highly  diverse  hydraulic  conditions  may  be  evaluated  as  poor  habitat.  Physical 
habitat  patterns  in  streams  can  be  described  more  realistically  in  a  methodo¬ 
logy  based  on  cumulative  frequencies  of  hydraulic  conditions.  Depiction  of  a 
stream  reach  in  terms  of  frequency  distributions  of  depth  and  velocity  is  more 
likely  to  capture  the  stream  heterogeneity  having  value  to  aquatic  biota  both  at 
the  population  and  community  levels  than  a  method  that  focuses  on  point-by- 
point  analyses.  If  a  cell  is  the  level  of  analysis  for  a  species,  it  seems 
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reasonable  that  the  appropriate  level  of  analysis  for  the  community  is  a 
frequency  analysis  of  cell  conditions  in  a  reach. 

The  holistic  perspective  of  the  RCHARC  provides  a  better  framework  to 
evaluate  the  ecosystem-level  differences  between  the  CSRS  and  the  project 
streams  by  better  describing  the  fluvial-geomorphic  factors  that  affect  fish 
habitat  (HUl,  Platts,  and  Beschta  1991).  The  system-level  perspective  is 
applied  later  in  this  report  using  the  Gavins  Point  tailwater  as  an  example  to 
better  detennine  how  the  project  alternatives  for  the  Missouri  River  differ  from 
the  CSRS.  The  stream  worker  can  use  this  approach  to  assess  broad  dif¬ 
ferences  in  hydraulic  patterns  between  the  CSRS  and  the  project  alternatives 
and  to  better  understand  how  patterns  of  physical  habitat  will  be  affected  by 
different  alternatives.  It  is  also  possible  to  gain  a  qualitative  understanding  of 
other  aquatic  processes  such  as  the  potential  for  bed  form  changes  or  potential 
for  organic  matter  transport. 

Rather  than  attempting  to  unweave  the  complex  tapestry  of  habitat  require¬ 
ments  for  each  species,  the  RCHARC  compares  the  underlying  patterns  of 
depths  and  velocities  in  the  two  systems  and  uses  the  results  as  the  basis  of  the 
community-level  impact  analysis.  The  degree  of  impact  is  approximated  by 
the  degree  to  which  the  physical  habitat  changes  between  the  target  and  stan¬ 
dard  systems.  Once  the  RCHARC  analysis  is  complete,  it  can  provide  infor¬ 
mation  useful  at  a  population  level  of  analysis  by  identifying  major  differences 
in  habitat  between  the  standard  and  target  systems.  These  differences  in  physi¬ 
cal  habitat  patterns  can  be  evaluated,  qualitatively  or  quantitatively,  to  deter¬ 
mine  habitat  features  whose  absence  may  account  for  the  dechne  of  endemic, 
rare,  or  endangered  species. 
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3  Implementing  the  RCHARC 


Steps  in  Applying  the  RCHARC 

The  RCHARC  is  composed  of  a  number  of  discrete  steps  summarized  in 
the  following  paragraphs.  Application  and  interpretation  of  the  RCHARC 
analysis  to  the  Gavins  Point  taUwater  are  discussed  in  the  succeeding  sections. 


Channel  descriptions 

The  CSRS  channel  and  project  alternative  channels  must  each  be  identified 
and  described  using  fine-resolution  cross-section  information.  Sufficient  infor¬ 
mation  for  the  CSRS  and  project  alternatives  must  be  available  to  allow 
description  of  discharge-specific  depth  and  velocity  distributions. 


Hydrologic  summaries 

The  dominant  hydrologic  patterns  associated  with  the  CSRS  and  each  of  the 
project  alternatives  must  be  summarized  using  summary  statistics  (e.g.,  median 
flows  or  flow  extremes)  and  time-steps  adequate  to  describe  the  hydrologic 
behavior  of  aU  systems. 


Stage-discharge  information 

Stage-discharge  information  for  the  CSRS  and  project  alternative  cross  sec¬ 
tions  must  be  obtained  from  hydraulic  models,  gauge  information,  or  field 
measurements.  The  stage-discharge  infonnation  is  combined  with  channel 
cross-section  information  to  generate  discharge-specific  depth  and  velocity 
distributions  for  the  CSRS  and  the  project  alternative  channels. 


Depth  and  velocity  distributions 

The  hydraulic  information  is  summarized  as  percent  distributions  of  depth 
and  velocity  for  both  the  CSRS  and  project  alternatives.  The  percent 
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distributions  for  depth  and  velocity  (treated  as  separate  variables  or  combined 
in  a  bivariate  depiction)  for  each  of  the  project  alternatives  are  evaluated  for 
their  similarity  to  the  CSRS  (Figure  3)  at  appropriate  time-steps  (e.g.,  monthly) 
using  qualitative  inspection,  correlation  analysis,  similarity  analysis,  or  other 
methods  for  determining  degree  of  similarity. 


Relative  impacts  of  each  project  alternative 

The  correlation  coefficients  for  each  project  alternative  are  summed,  month 
by  month,  over  the  period  of  the  hydrologic  record.  The  project  alternative 
that  provides  depth  and  velocity  conditions  most  nearly  like  the  CSRS  will 
have  the  largest  sum,  and  the  project  alternative  having  the  least  similarity 
should  have  the  smallest  sum. 


Correct  for  top  width 

In  some  cases  it  may  be  possible  for  the  top  widths  of  the  CSRS  and  proj¬ 
ect  alternative  chaimels  to  vary  substantially  even  though  the  depth  or  velocity 
distributions  are  very  similar.  As  an  option,  the  coirelation  coefficients  can  be 
linearly  adjusted  by  the  degree  to  which  project  alternative  channel  top  widths 
differ  from  the  CSRS  channel  topwidths. 


Trends  and  patterns  in  the  results 

As  with  any  methodology  that  distiUs  complex  systems  into  single  number 
rankings,  it  is  important  to  explore  the  significance  of  differences  and  similari¬ 
ties  among  project  alternatives  and  between  the  project  alternatives  and  the 
CSRS. 


Channel  Descriptions 

Physical  habitat  in  the  CSRS  and  the  project  alternatives  must  be  described 
using  channel  cross-section  data.  Transect  data  must  be  of  sufficient  detail  to 
allow  adequate  characterization  of  aquatic  habitat — usually  requiring  at  least 
20  points  in  the  channel  and  a  vertical  accuracy  of  0.5  ft'  for  a  system  as 
large  as  the  Missouri  River.  A  sufficient  numter  of  transects  must  be  identi¬ 
fied  and  selected  to  characterize  the  habitat  variability  of  a  river  reach. 
Usually,  a  minimum  of  four  to  eight  transects  is  required. 

MRD  has  collected  and  archived  channel  cross-section  data  for  the 
Missouri  River  at  approximately  5-year  intervals  beginning  with  cross 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  viii. 
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examples  would  have  more  pairs).  The  H1/Q1  pair  in  conjunction  with  other  information  is  used  to  determine  the  distribution  of 
depth  D  or  velocity  V  in  the  channel.  A  similar  process  is  followed  to  determine  the  depth  or  velocity  distribution  in  the  project 
alternative.  The  resulting  distributions  are  compared  using  correlation  analysis  or  a  similarity  statistic  to  determine  how  a 
project  channel  compares  to  a  standard  channel.  This  process  would  be  repeated  nine  times  to  generate  the  table  of 
correlation  coefficients  in  the  figure 


sections  collected  prior  to  or  immediately  after  closure  of  each  of  the  dams  and 
extending  to  the  present.  Each  cross  section  is  composed  of  approximately 
100  points  with  a  vertical  accuracy  of  0.1  to  0.2  ft.  Available  channel  cross 
sections  were  evaluated  by  MRD.  Between  6  and  10  cross  sections  for  each 
tailwater  were  selected  to  represent  major  habitats  in  the  river  from  the  dam  to 
the  headwaters  of  the  next  downstream  reservoir.  Major  habitats  included 
backwaters,  chutes,  main  channels,  and  sandbar  areas.  For  each  tailwater, 
channel  cross  sections  included  major  habitat  features  in  approximately  the 
same  proportion  that  they  existed  in  the  entire  river  reach.  State  and  Federal 
biologists  used  aerial  photographs,  existing  cross  sections,  and  site  visits  to 
select  transects  for  measurement.  Preregulation  transects  were  selected  from 
historic  data  collected  at  the  same  locations  or  as  close  as  possible  to  the  trans¬ 
ects  selected  for  measurement. 

Cross  sections  for  Gavins  Point  were  categorized  by  top  width  as  narrow, 
wide,  transitional,  or  divided.  After  selection  of  cross  sections  from  MRD 
archives,  each  location  was  resurveyed  under  both  a  high-  and  low-flow  condi¬ 
tion  using  standard  stream  gauging  methods  to  define  the  channel  conditions 
for  project  alternatives.  Results  of  the  two  surveys  were  used  to  prepare  two 
velocity  calibration  data  sets  that  would  be  used  to  calibrate  the  hydraulic 
models  used  to  simulate  lateral  cross-section  velocities.  Two  resurveys  were 
used  to  get  a  better  approximation  of  the  channel  shape  associated  with  the 
two  dominant  flows.  Median  flow  used  to  support  navigation  during  the  navi¬ 
gation  season  was  approximately  30,000  cfs.  Median  flow  during  the  non¬ 
navigation  season  was  approximately  10,000  cfs. 

Channel  cross-section  and  velocity  data  for  this  predominantly  sand  bed 
river  were  also  collected  at  a  high  flow  and  a  low  flow  because  the  channel 
shape  might  change  between  the  two  flow  extremes.  This  reduced  the  risk  that 
a  channel  cross  section  obtained  under  high-flow  conditions  would  not  be 
representative  of  channel  shape  at  lower  discharges,  particularly  at  the  edges  of 
the  channel.  However,  comparison  of  the  high-  and  low-flow  cross  sections  at 
both  flows  at  similar  locations  indicated  only  relatively  minor  changes  in  chan¬ 
nel  shape.  The  two  velocity  calibration  data  sets  were  collected  near  the  two 
end  points  of  the  flow  range  representative  of  most  of  the  operational  alterna¬ 
tives,  ensuring  that  any  change  in  flow  pattern  would  be  captured  by  the 
two  calibration  data  sets. 

The  complete  channel  data  set  available  for  RCHARC  analysis  consisted  of 
river  cross-section  information  collected  at  the  same  locations  (or  usually 
within  several  tenths  of  a  mile)  for  Gavins  Point  tailwater  prior  to  or  near 
closure  of  the  dam  immediately  upstream  and  again  in  1989  or  more  recently. 
Detailed  descriptions  of  the  channel  cross-section  information  available  are 
described  later  in  this  section. 

The  Gavins  Point  cross  sections  were  resurveyed  13-19  March  and 
10-19  July  1989  by  a  team  composed  of  representatives  of  the  Nebraska  Game 
and  Parks  Commission,  the  U.S.  Geological  Survey  (USGS),  and  U.S.  Army 
Corps  of  Engineers.  Flows  ranged  from  8,000  to  10,500  cfs  during  the  March 
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data  collection  and  were  approximately  30,000  cfs  during  the  July  measure¬ 
ments.  Eight  transects  in  the  reach  were  selected  for  data  collection,  each 
representing  a  different  river  type.  Only  six  locations  were  surveyed  during 
March  because  of  bad  weather  and  increased  releases  from  Gavins  Point  for 
the  start  of  the  navigation  season.  In  March,  cross-section  data  were  collected 
at  RM  778.90,  780.92,  783.61,  786.73,  797.50,  and  804.28.  July  cross  sections 
were  at  755.56,  778.90,  780.92,  783.61,  786.73,  793.60,  797.50,  and  804.28. 
The  preproject  cross  sections  were  collected  from  1955  through  1959. 


Hydrologic  Summaries 

The  description  of  the  GSRS  consists  of  both  channel  characteristics  occur¬ 
ring  over  long-term  discharge  patterns.  River  discharges  vaiy  by  location  on 
the  Missouri  River  because  of  tributary  inflow.  Monthly  discharges  for  the 
period  of  record  to  characterize  the  GSRS  were  obtained  from  LRS  model 
output  using  a  run  in  which  the  storage  of  the  reservoirs  was  “held  out” 
Therefore,  the  ou^ut  of  the  LRS  model  would  simulate  the  hydrologic  behav¬ 
ior  of  the  Missouri  River  as  though  the  dams  had  not  been  constructed.  These 
“holdout  runs”  were  used  to  generate  preproject  hydrological  patterns. 

Use  of  the  LRS  model  allowed  MRD  to  simulate  the  preproject  hydrology 
of  the  Missouri  River  before  data  were  consistently  available  for  all  present 
gauge  locations;  that  is,  the  LRS  model  was  able  to  synthesize  gauge  informa¬ 
tion  at  locations  having  incomplete  records.  Additionally,  the  LRS  model 
could  simulate  the  period  of  record  after  the  project  was  constructed  as  though 
the  dam  did  not  exist.  Use  of  existing  gauge  information  to  characterize  the 
preproject  Missouri  River  provided  a  hydrologic  record  limited  to  a  duration  as 
short  as  24  years  for  one  of  the  main  stem  Missouri  River  tailwaters.  A 
24-year  record  is  inadequate  to  describe  relatively  long-term  drought-flood 
cycles  that  characterize  the  Missouri  River.  The  LRS  model  could  provide  a 
93-year  hydroperiod  to  more  completely  characterize  the  preproject 
hydrological  conditions. 

Variation  in  flows  from  year  to  year  can  have  important  influences  on  the 
distribution  of  aquatic  organisms.  Variation  in  flows  may  favor  species  that 
evolved  in  the  uruegulated  river  system  under  naturally  occurring  discharge 
and  channel  conditions  at  the  expense  of  exotic  or  introduced  species.  Gonse- 
quently,  consideration  of  the  median  or  mean  flow  only  may  not  adequately 
relate  the  effects  of  the  preproject  hydrology  on  the  integrity  of  the  native 
warmwater  fish  community.  Three  categories  of  preproject  flows  were  each 
considered  to  define  a  GSRS;  low  flow  (75  percent  exceedance),  median  flow 
(50  percent  exceedance),  and  high  flow  (25  percent  exceedance)  conditions  for 
the  Gavins  Point  reach,  as  sho\vn  in  the  following  tabulation: 
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Month 

Median 

High 

Low 

1 

12.0 

15.0 

9.0 

2 

16.5 

21.0 

13.5 

3 

40.0 

50.0 

32.0 

4 

38.0 

49.0 

27.0 

5 

47.0 

59.0 

34.0 

6 

81.0 

106.0 

55.0 

7 

42.0 

59.0 

29.0 

8 

19.0 

24.0 

14.0 

9 

17.0 

23.0 

3.5 

10 

17.5 

22.0 

14.0 

11 

16.5 

20.0 

13.5 

12 

11.0 

14.0 

8.0 

The  annual  monthly  time  series  associated  with  each  flow  category  for  each 
year  was  derived  firom  analysis  of  data  generated  by  the  LRS  model  and  pro¬ 
vided  to  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES)  by 
MRD.  An  artificial  “median  flow  year”  was  composed  of  the  median  flow 
from  each  month.  For  example,  the  median  flow  year  for  the  Gavins  Point 
tailwater  was  synthesized  by  obtaining  the  50  percent  exceedance  flow  for 
each  month  for  the  Yankton  Gage  (a  gage  in  the  tailwater  of  Gavins  Point).  A 
conceptually  similar  process  was  used  to  synthesize  the  low-flow  and  high- 
flow  CSRS’s  for  the  Gavins  Point  tailwater.  During  application  of  the 
RCHARC  to  the  93-year  period  of  hydrologic  record,  each  water  year  under 
the  holdout  operation  of  the  LRS  was  assigned  to  one  of  the  three  exceed¬ 
ance  categories  according  to  rules  developed  by  the  MRD. 


Hydraulic  Information 

Lateral  flow  pattern  descriptions 

IFG-4,  one  of  the  hydraulic  programs  within  the  PHABSIM  system,  was 
used  to  develop  the  lateral  flow  (velocity)  pattern  at  each  transect  using  infor¬ 
mation  from  stage-discharge  relationships,  channel  cross-section  information, 
and  two  velocity  calibration  data  sets.  The  hydraulic  component  of  the 
PHABSIM  system  assumes  that  the  shape  of  the  channel  does  not  change  with 
streamflow  over  the  range  of  flows  being  simulated.  The  results  of  the 
hydraulic  calculations  are  water  surface  elevations  and  velocities.  The  water 
surface  elevations  are  one-dimensional  in  that  the  same  value  is  used  for  any 
point  on  a  cross  section.  In  contrast,  the  velocity  varies  from  point  to  point 
across  any  cross  section. 
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Several  options  within  the  IFG-4  program  are  available  to  generate  cell 
depths  and  velocities  for  a  given  discharge.  The  IFG-4  program  can  employ 
several  methods  to  synthesize  cell-by-ceU  velocities  if  cdibration  velocities  are 
unavailable.  For  this  application,  the  hydraulic  radius  option  was  used  to 
separate  the  total  discharge  at  each  transect  into  cells.  This  option  results  in 
deeper  cells  having  greater  water  velocities.  Alternatively,  velocity  measure¬ 
ments  made  in  the  field  under  steady  flows  at  one  or  more  discharges  can  be 
used  for  calibration.  The  measured  velocities  are  used  to  solve  Manning’s 
equation  for  Manning’s  n  for  each  cell.  The  calculated  cell-specific  n  values 
are  then  used  to  generate  estimates  of  velocities  for  each  cell  over  a  range  of 
simulated  discharges.  After  estimating  a  lateral  velocity  pattern  based  on  the 
calculated  Manning’s  n,  the  IFG-4  program  then  checks  Ae  calculated  water 
surface  elevation  against  the  given  water  surface  elevation  provided  from 
another  source  (described  in  the  stage-discharge  section  above)  and  then,  if 
necessary,  modifies  all  cell  velocities  by  a  common  factor  to  raise  or  lower  the 
estimated  water  surface  elevation  until  it  matches  the  given  water  surface 
elevation. 

The  flow  pattern  often  changes  as  channel  discharge  changes.  Conse¬ 
quently,  it  is  important  to  select  either  an  intermediate  discharge  for  the  veloc¬ 
ity  calibration  data  set  that  can  be  applied  to  a  relatively  narrow  range  of  flows 
or  to  select  two  or  mote  velocity  calibration  data  sets  depending  upon  the 
discharge  range  being  simulated  and  the  rate  at  which  the  flow  pattern 
changes.  For  the  Missouri  River  application,  two  velocity  calibration  data  sets 
were  collected,  one  near  the  upper  limit  of  normal  operation  (navigation  season 
releases — April  through  November)  and  one  near  the  lowest  discharge  (non¬ 
navigation  season  releases — ^December  through  March)  that  would  be  ordi¬ 
narily  released  under  most  alternatives.  Calibration  of  the  IFG-4  program  to 
measured  velocities  is  generally  superior  to  using  hydraulic  radius  for  estimat¬ 
ing  lateral  flow  patterns  (Milhous,  Updike,  and  Schneider  1989). 

The  lower  velocity  calibration  discharge  was  applied  from  the  lowest  simu¬ 
lated  stage-discharge  pair  to  about  the  midpoint  of  the  flow  difference  between 
the  two  calibration  data  sets.  For  the  Gavins  Point  application,  the  low-flow 
calibration  data  set  collected  at  about  11,000  cfs  was  applied  from  6,000  to 
18,(X)0  cfs.  The  high-flow  velocity  calibration  data  set  was  applied  from  the 
midpoint  discharge  to  the  highest  discharge  simulated.  For  the  Gavins  Point 
application,  the  high-flow  velocity  calibration  data  set  collected  at  about 
30,000  cfs  was  applied  from  18,000  to  50,000  cfs.  The  preproject  channel 
cross  sections  did  not  include  any  velocity  measurements.  Channel  velocities 
were  predicted  by  the  IFG-4  program  based  on  hydraulic  radius  for  the  pre¬ 
project  channel  cross  sections.  ' 

Channel  geometry,  velocity,  substrate,  and  discharge  information  were 
provided  by  the  MRD  from  field  surveys.  From  this  information,  two  separate 
IFG-4  input  data  sets  were  developed.  One  IFG-4  input  data  set  was  devel¬ 
oped  for  the  discharges  covered  by  the  low-flow  velocity  calibration  data  sets 
and  one  set  for  the  discharges  covered  by  the  high-flow  velocity  calibration 
data  sets.  The  IFG-4  program  generates  a  binary  output  data  set  for  each  input 
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data  set.  The  IFG-4  output  data  sets  are  composed  of  tables  of  cell-by-cell 
predicted  depths  and  velocities  for  each  cross  section  and  discharge.  The 
IFG-4  binary  ouq)ut  data  sets  were  converted  to  ASCII  using  the  LSTVDX 
program  of  the  PHABSIM  system.  The  two  output  data  sets  associated  with 
each  transect  were  appended  after  processing  by  the  LSTVDX  program. 
Standard  options  were  used  in  the  IFG-4  program  except  that  the  HABTAM- 
HABTAV  method  of  using  vertical  velocities  instead  of  average  cell  velocities 
was  employed.  Use  of  average  cell  velocities  has  a  tendency  to  reduce  the 
variability  in  the  data  set  by  reducing  high  velocities  or  increasing  lower  veloc¬ 
ities  by  averaging  with  less  extreme  neighboring  cells.  IFG-4  was  run  for  the 
sequences  of  discharges  listed  in  the  following  tabulation: 


Discharge,  cfs 

Preregulation 

High 

Low 

6,000 

20,000 

6,000 

8,000 

24,000 

8,000 

10,000 

28,000 

10,000 

12,000 

32,000 

12,000 

14,000 

36,000 

14,000 

16,000 

40,000 

16,000 

20,000 

46,000 

24,000 

50,000 

— 

28,000 

32,000 

36,000 

40,000 

46,000 

50,000 

An  example  of  IFG-4  output  processed  by  the  LSTVDX  program  is  pre¬ 
sented  in  Table  1.  Three  different  sets  of  IFG-4  data  sets  were  processed. 
Preproject  IFG-4  data  sets  were  generated  in  which  the  velocities  were  based 
on  hydraulic  radius  since  measured  velocities  were  unavailable  for  these  trans¬ 
ects.  Two  postproject  data  sets  were  created.  The  first  postproject  data  set 
had  the  calibration  velocities  stripped  out,  and  the  IFG-4  program  generated 
output  velocities  based  on  hydraulic  radius.  The  second  postproject  data  set 
included  calibration  velocities,  and  the  IFG-4  program  predicted  output  veloci¬ 
ties  based  on  measured  velocities.  For  the  second  postproject  data  set,  veloc¬ 
ities  were  based  on  cell-specific  Manning’s  n  obtained  from  the  two  velocity 
calibration  data  sets. 
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Depth  and  velocity  distributions 

The  three  ASCII-format  tables  of  depth  and  velocity  generated  for  each 
transect  were  then  evaluated  statistically  using  the  Statistical  Analysis  System 
(SAS  Institute,  Inc.,  1988).  The  following  steps  were  involved  in  the  statistical 
analysis: 

a.  Cell-by-cell  information  obtained  from  the  LSTVDX  program  was  pro¬ 
cessed  to  eliminate  cells  having  zero  depths  (cells  above  the  water 
surface)  for  each  transect. 

b.  Cells  from  similar  cross-section  categories  were  combined;  e.g.,  all  cells 
from  wide  transects  were  placed  into  a  single  category  and  individual 
cross  sections  were  no  longer  distinguished. 

c.  Depths  were  rounded  to  the  nearest  0.5  ft  and  velocities  were  rounded 
to  the  nearest  0.25  fps  using  the  ROUND  function  of  SAS  (SAS 
Institute,  Inc.,  1988). 

PROC  FREQ  of  the  SAS  was  employed  to  generate  percent  cumulative 
frequencies  of  depth  and  velocity  for  each  of  the  three  IFG-4  output  data  sets 
for  each  discharge  and  each  separate  category  of  transects.  Separate  frequen¬ 
cies  were  generated  for  depth  and  velocity  for  the  univariate  application.  For 
the  bivariate  application,  a  single  frequency  distribution  was  generated  that 
combined  both  variables. 


Univariate  Analysis 

Depth-velocity  distribution 

Inspection  of  the  cumulative  frequency  distributions  generated  by  the  fre¬ 
quency  analysis  indicated  that  they  closely  followed  the  general  pattern  exhib¬ 
ited  by  the  logistic  equation.  Frequency  distributions  for  depth  and  velocity 
for  each  of  the  three  IFG-4  output  data  sets  (preproject,  project  without  calibra¬ 
tion  velocities,  and  project  with  calibration  velocities)  for  discharge  and  chan¬ 
nel  category  were  then  fit  to  the  logistic  equation 

^  (1) 


1  + 


K-F 


where 

CPD  =  cumulative  percent  composition 

K  =  apparent  maximum  cumulative  percent  (usually  near  100) 
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F  =  percent  composition  of  the  first  depth  or  velocity  increment 

R  =  rate  of  change  from  one  cumulative  depthA'clocity  increment  to  the 
next 

X  =  midpoint  of  deptlVvelocity  increment 

The  fit  of  the  logistic  equation  to  the  cumulative  frequency  distributions 
(one  distribution  for  each  discharge  and  channel  category)  was  highly  signifi¬ 
cant.  Table  2  presents  a  summary  of  the  statistics  for  the  fit  of  the  logistics 
equation  for  the  Gavins  Point  tailwater. 

Curve-fitting  was  employed  for  three  reasons.  First,  the  lower  Missouri 
River  is  a  sand  bed  stream  characterized  by  zones  of  uniform  depths.  How¬ 
ever,  field  observation  indicates  that  these  uniform  depths  may  gradually 
increase  or  decrease  in  either  the  upstream  or  downstream  direction.  Distribu¬ 
tions  of  depth  and  velocity  based  on  curve-fitted  data  would  tend  to  smooth 
out  the  effect  of  a  particularly  abundant  depth  or  velocity  category.  These 
distributions  would  thus  be  less  likely  to  be  influenced  by  the  localized  charac¬ 
teristics  of  a  specific  location  on  the  river  and  be  more  likely  to  represent  a 
generalized  feature  of  the  river.  This  was  a  particular  concern  for  the  Gavins 
Point  reach  data,  which  had  relatively  few  cross  sections  and  could  be  easily 
biased  by  a  high  percentage  of  a  particular  depth  or  velocity  category  in  the 
limited  selection  of  cross  sections. 

Second,  curve-fitting  allowed  the  use  of  functions  to  describe  each  of  the 
distributions.  Functions  can  be  stored,  adjusted,  and  manipulated  more  effi¬ 
ciently  than  raw  data.  Third,  with  curve-fitting,  the  cumulative  frequency 
distribution  of  intennediate  discharges  (discharge  increments  not  simulated) 
could  be  estimated  by  linear  interpolation  of  the  coefficients  between  neighbor¬ 
ing  discharges. 

Corrections  to  the  preproject  velocity  predictions  were  made  by  estimating 
the  error  associated  with  use  of  hydraulic  radius  to  predict  velocity  distribu¬ 
tions  on  project  data.  The  project  IFG-4  data  sets  can  be  used  to  evaluate  the 
error  associated  with  the  use  of  hydraulic  radius  to  estimate  velocity  by  mn- 
ning  the  IFG-4  program  with  and  without  the  velocity  calibration  data  sets. 

The  differences  (residuals)  between  the  velocity  cumulative  frequency  distribu¬ 
tions  represent  the  error  in  use  of  the  hydraulic  radius  to  predict  velocity. 
Inspection  of  the  patterns  in  the  residuals  indicated  that  each  could  be  fit  by  a 
quadratic  equation. 

It  was  decided  that  adding  back  the  residuals  in  the  form  of  a  quadratic 
equation  (Figures  4-6)  to  the  preproject  data  set  would  provide  the  most  accu¬ 
rate  estimate  of  the  preproject  cumulative  velocity  distribution  because  the  pre- 
project  cross  sections  were  at  or  very  near  the  locations  of  the  project  cross 
sections,  the  river  had  not  moved  or  meandered  substantially,  and  channel  top 
widths  had  not  changed  dramatically. 
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Figure  4.  Depiction  of  method  used  to  correct  the  CSRS  velocity  distributions.  The  cumula¬ 
tive  frequency  distribution,  represented  by  the  logistic  equation,  is  corrected  using  a 
quadratic  equation  derived  from  analysis  of  project  data  with  and  without  velocity 
calibration  data.  A  frequency  distribution  is  then  obtained  from  the  corrected  cumu¬ 
lative  frequency  distribution 


If  the  study  objective  had  been  restricted  to  the  ranking  of  alternatives  to  a 
standard,  velocities  based  on  hydraulic  radius  would  have  been  used  for  both 
the  CSRS  and  project  channels.  The  same  velocity  bias  would  occur  in  both 
the  CSRS  and  the  project  alternatives  and,  therefore,  the  relative  ranking  of 
alternatives  should  not  be  affected.  However,  the  results  from  this  study  may 
also  be  employed  for  river  restoration,  and  thus  the  velocities  predictions 
should  be  as  accurate  as  possible. 

At  this  point  in  the  analysis  using  the  RCHARC,  a  discharge-specific  and 
channel  category-specific  cumulative  frequency  distribution  for  depth  and 
velocity  (corrected  as  described  above)  for  the  preproject  transects  was  avail¬ 
able  along  with  corresponding  information  for  the  project  channel  alternatives. 
Simple  frequency  distributions  were  derived  from  die  cumulative  frequency 
distributions  by  subtracting  from  each  cumulative  depth  or  cumulative  velocity 
increment  the  cumulative  increment  immediately  preceding  it  (e.g..  Figures  7 
and  8).  The  final  product  of  this  step  was  discharge-specific  simple  frequency 
distributions  of  depth  and  velocity  for  each  channel  category  within  the  tail- 
water.  Additionally,  by  summing  the  top  widths  of  each  cell  for  each  channel 
category  within  the  tailwater,  it  was  also  possible  to  estimate  the  discharge- 
specific  top  widths  associated  with  the  preproject  and  project  channels.  Dis¬ 
charges  were  rounded  to  the  nearest  thousand  cubic  feet  per  second. 


Relative  impacts  of  each  project  alternative 

The  depth  and  velocity  distributions  associated  with  each  project  alternative 
were  expressed  differently  from  the  depth  and  velocity  patterns  associated  with 
the  CSRS’s.  Depth  and  velocity  distributions  associated  with  each  alternative 
were  not  assessed  directly,  but  rather  interpolated  from  incremental  sequences 
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Figure  6.  Comparison  of  velocity  distributions  made  at  a  50,000-cts  discharge  (symbols  as  defined  for  Figure  5) 


Figure  7.  Comparison  of  the  depth  and  velocity  frequencies  between  the  CSRS  at  a  dis¬ 
charge  of  18,000  cfs  (historical)  and  project  channel  conditions  at  a  discharge  of 
28,000  cfs  (operational)  for  the  wide  channel  category.  Note  the  greater  depths 
and  lower  velocities  in  the  historical  channel  than  in  the  project  channel.  The 
velocity  distributions  and  depth  distributions  do  not  necessarily  have  to  mass 
balance  if  the  top  widths  are  substantially  different 
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Figure  8.  Comparison  of  the  depth  and  velocity  frequencies  between  the  CSRS  at  a  dis¬ 
charge  of  43,000  cfs,  project  channel  at  a  discharge  of  40,000  cfs,  and  special 
holdout  run  to  simulate  a  year  as  though  the  dams  were  not  in  place 


of  nows  for  each  site.  For  example,  project  alternative  Hows  were  simulated 
for  6,000,  8,000,  10,000,  12,000,  14,000,  16,000,  20,000,  24,000,  28,000, 
32,000,  36,000,  40,000,  46,000,  and  50,000  cfs  for  the  Gavins  Point  reach. 

The  lower  How  limit  was  determined  by  the  limitations  of  the  HEC-2  model  to 
simulate  Hows  less  than  6,000  cfs,  and  the  upper  limit  was  determined  by  the 
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discharge  at  which  overbank  flow  occurred.  The  MRD  felt  that  stages  for 
overbank  flows  could  not  be  accurately  simulated.  The  CSRS’s  were 
characterized  by  distributions  of  depth  and  velocity  associated  with  specific 
water  years  on  a  month-by-month  basis.  For  example,  the  depth  distribution 
for  January  (1 1,000  cfs)  in  a  median  water  year  was  obtained  by  linear  inter¬ 
polation  of  the  and  A"  coefficients  of  the  logistic  equation  for  10,000  and 
12,000  cfs.  The  water- year-specific  CSRS  was  described  in  terms  of  the  depth 
and  velocity  distribution  for  each  month. 

The  relative  value  of  each  alternative  flow  at  a  particular  month  for  one  of 
the  three  CSRS’s  was  obtained  by  detemiining  the  similarity  between  the 
depth/velocity  patterns  associated  with  monthly  CSRS  flows  and  the 
depth/velocity  increments  used  to  evaluate  the  postproject  channel  (detailed 
description  of  the  steps  is  shown  in  Figures  4,  9,  10,  and  1 1).  The  impact  of  a 
particular  operational  alternative  for  1  year  could  be  determined  by  looking  up 
its  similarity  to  preproject  distributions  on  a  table  of  incremental  project  dis¬ 
charges  versus  preproject  discharges  (Appendix  A).  Each  of  the  12  resulting 
similarities  is  summed  to  determine  the  impact  of  that  alternative.  The  process 
can  be  expanded  to  include  any  number  of  months.  The  correlation  coeffi¬ 
cients  were  rescaled  to  range  from  0.0  to  1.0  instead  of  -1.0  to  1.0. 

A  variety  of  similarity  metrics  can  be  employed;  however,  for  this  analy¬ 
sis,  the  Pearson  product-moment  correlation  analysis  was  used  to  relate  how 
the  depth  or  velocity  categories  varied  between  the  CSRS  and  a  series  of  incre¬ 
mental  project  flows.  Other  statistics  were  evaluated  at  the  beginning  of  this 
study,  but  the  ranking  of  coefficients  of  flow  increments  versus  a  particular 
historic  flow  did  not  change.  Also,  as  pointed  out  by  Clifford  and  Stephenson 
(1975),  most  methods  of  determining  similarity  produce  the  same  general 
patterns  of  results.  However,  further  use  of  the  RCHARC  may  determine  that 
other  statistics  are  more  suitable. 


Correction  for  top  width 

It  is  possible  to  have  similar  depth  or  velocity  distributions  (habitat  quality) 
but  very  different  top  widths  (habitat  quantity).  To  account  for  differences  in 
quantity  between  the  CSRS  and  the  project  alternatives,  the  correlations 
between  the  CSRS’s  and  the  project  alternatives  were  adjusted  by  a  coefficient 
keyed  to  the  difference  between  the  preproject  top  width  and  the  top  widths 
associated  with  each  alternative  flow  increment.  For  example,  if  the  top  width 
of  a  historical  flow  was  identical  to  a  particular  flow  increment,  then  there 
would  be  no  adjustment  of  the  correlation  coefficient.  However,  if  the  pre¬ 
project  top  width  at  a  flow  associated  with  a  particular  month  was  less  than  a 
particular  flow  increment,  then  the  correlation  coefficient  was  reduced  by  the 
ratio  difference  between  the  two  top  widths  as 
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discharges.  The  cumulative  frequency  is  fitted  to  the  logistic  equation  and  a  smoothed  velocity  and  depth  distribution  is 
developed  from  each  fitted  cumulative  frequency  distribution.  In  this  case,  the  analysis  is  performed  for  only  1  month 
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weighted  more  heavily  if  necessary 


SUM  WEIGHTED 
CORRELATION  COEFFICIENTS 


HISTORICAL 


HISTORICAL 


ALTA 

+  0.04 

+  0.17 

ALT  B 

+  0.04 

+  0.16 

ALT  C 

+  0.08 

+  0.16 

HISTORICAL 

DIVIDED 


HISTORICAL 


NARROW 


HISTORICAL 


ALTA 

+  0.15 

+  0.19 

ALT  B 

+  0.15 

+  0.21 

ALTC 

+  0.07 

+  0.21 

TRANSITIONAL 


WIDE 


ALTA 
ALT  B 
ALT  C 


+  0.55 


+  0.56 


+  0.52 


Figure  1 1 .  Each  of  the  weighted  correlation  coefficients  from  the  rightmost  column  of  Figure  10 
is  then  summed  to  represent  a  single  value  of  impact  for  each  of  the  three  alterna¬ 
tives.  For  this  example,  covering  Figures  4,  9,  10,  and  this  figure,  the  example  was 
performed  for  only  1  month.  For  multiple  months  (or  other  time-steps),  the  correla¬ 
tion  coefficients  would  be  summed  across  rows  to  develop  a  single  value  of  impact 


COR  COEF  =  COEFA* 


_  PROJJTOP  -  PREJTOP^ 

PROJ  TOP 

—  J 


(2) 


where 

COR  COEF  =  adjusted  correlation  coefficient 
COEFA  =  unadjusted  correlation  coefficient 
PRE  TOP  =  CSRS  (preproject)  top  width 
PROJJTOP  =  project  top  width 

If  the  preproject  top  width  at  a  flow  associated  with  a  particular  month  was 
greater  than  a  particular  flow  increment,  then  the  correlation  coefficient  was 
computed  as 


f 

COR_COEF  =  COEFA*  1 

V 


■x 

PREjrOP  -  PROJJTOP 

PRE  TOP 

—  J 


(3) 
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For  Gavins  Point  tailwater,  48  tables  of  correlation  matrices  were  generated, 
categorized  as  follows: 

a.  Three  different  water  years  (median,  low,  and  high). 

b.  Two  flow  variables  (depth  and  velocity). 

c.  Four  channel  categories  (divided,  narrow,  transitional,  and  wide). 

d.  Two  top  width  categories  (adjusted  and  unadjusted). 

The  different  channel  categories  for  Gavins  Point  could  be  combined  by 
weighting  each  category  by  the  percent  distance  of  the  reach  that  it  represented 
and  summing  the  correlation  coefficients  (Figure  11). 


Bivariate  Analysis 

Depth-velocity  distributions 

The  following  steps  were  involved  in  creating  bivariate  data  sets  of  depth- 
velocity  profiles  for  the  Gavins  Point  reach  of  the  Missouri  River; 

a.  Cell-by-cell  information  obtained  from  the  LSTVDX  program  was  pro¬ 
cessed  to  eliminate  cells  above  the  water  surface  for  each  transect. 

b.  Cells  from  similar  cross-section  categories  were  combined  (as  in  the 
univariate  analysis),  and  individual  cross  sections  were  no  longer  distin¬ 
guished.  In  addition  to  the  channel  categories  (divided,  narrow,  transi¬ 
tional,  and  wide),  transects  were  classified  as  to  proximity  to  the  dam 
(closer  than  20  miles  is  near;  farther  than  20  miles  is  far).  Also,  chan¬ 
nel  categories  were  combined  for  analyzing  overall  trends  within  the 
Gavins  Point  reach. 

c.  Depth  and  velocity  values  were  converted  to  metric  (1  ft  =  0.3048  m) 
and  rounded  off  (depth  to  nearest  0.5  m,  velocity  to  nearest  0.2  m/sec). 
Discharge  alternatives  were  rounded  to  the  nearest  1,000  cfs. 

Bivariate  frequency  distributions  (depth-velocity)  per  discharge  were  gener¬ 
ated  using  the  SAS  FREQ  procedure,  with  each  depth-velocity  increment 
across  each  section  weighted  by  the  water  surface  area  (cell  width)  it  repre¬ 
sented.  The  sparse  option  was  used  in  the  FREQ  procedure  to  write  out  all 
possible  depth- velocity  combinations.  Frequency  distributions  were  created  in 
this  manner  for  aU  discharges  per  channel  category  for  both  project  and  pre¬ 
project  data  sets. 

Rounding  the  depth  and  velocity  values  resulted  in  a  spiking  of  the  fre¬ 
quency  distribution  data  at  the  rounded  increments.  To  alleviate  this  inherent 
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rounding  problem,  the  sharp  peaks  in  the  bivariate  data  set  were  smoothed  by 
taking  a  running  average  of  three  values  of  the  percent  frequency  variable. 


Relative  Impacts  of  each  project  alternative 

Project  alternative  flows  were  simulated  for  the  same  14  discharges  as  in 
the  univariate  analysis,  with  historical  depth-velocity  frequency  distributions 
interpolated  from  incremental  sequences  of  flows.  The  CSRS’s  were  charac¬ 
terized  by  bivariate  distributions  associated  with  monthly  flows  for  a  median 
water  year. 

To  examine  the  extent  of  similarity  (or  dissimilarity)  between  project  and 
preproject  depth-velocity  profiles,  the  quantitative  Canberra  metric  coefficient 
(Lance  and  Williams  1966,  1967)  was  employed.  The  Canberra  coefficient 
technique  was  chosen  because  its  resemblance  value  reflects  high  and  low 
scoring  attributes  evenly.  The  metric  coefficient  is  based  on  the  average  of  a 
series  of  fractions  relating  to  interentity  resemblance,  thus  having  inherent 
attribute  standardization  (Boesch  1977).  The  dissimilarity  Canberra  metric 
coefficient  has  the  form: 


-OWE  K; 


(4) 


where 

Aji^  =  Canberra  metric  coefficient 

n  =  total  number  of  attributes  excluding  aU  double  zero  matches 

X-j  =  preproject  frequency  distribution  per  depth-velocity  increment 

Xj^  =  project  frequency  distribution  per  depth-velocity  increment 

In  cases  where  one  attribute  within  a  match  is  zero,  a  value  of  0.01  is  sub¬ 
stituted  to  ensure  that  varying  attribute  differences  are  reflected  accordingly. 
The  Canberra  coefficient  is  defined  on  the  range  0.0  <  <  1.0,  where  = 

1.0  indicates  maximum  dissimilarity.  To  convert  the  metric  to  a  similarity 
measure,  simply  subtract  Ay^  from  1.  All  Canberra  coefficients  reported  here 
have  been  converted  to  measures  of  similarity  to  match  the  correlation  coeffi¬ 
cients  of  the  univariate  analysis. 

As  in  the  univariate  analysis,  the  difference  in  top  width  between  the  histor¬ 
ical  conditions  and  the  project  alternatives  was  accounted  for  when  describing 
the  degree  of  similarity  between  the  two  data  sets.  Refer  to  Part  III  for  the  top 
width  adjustment  equation. 
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For  the  Gavins  Point  reach,  Canberra  similarity  matrices  were  generated  in 
the  bivariate  analysis  for  a  median  water  year,  categorized  as  follows: 

a.  Four  channel  categories  (divided,  narrow,  transitional,  and  wide). 

b.  Two  top  width  categories  (adjusted  and  unadjusted). 

To  determine  impact  for  the  Gavins  Point  reach  at  a  single  discharge,  the 
four  channel  categories  were  combined  by  weighting  each  category  by  the 
proportional  distance  of  the  tailwater  that  it  represented  and  summing  the  Can¬ 
berra  coefficients  (Table  3).  Monthly  median  historical  deptlWelocity  fre¬ 
quency  distributions  were  linearly  interpolated;  i.e.,  median  discharge  for 
April  =  38,000  cfs,  so  frequency  distributions  were  calculated  using  discharges 
of  36,000  and  40,000  cfs.  Top  widths  associated  with  median  monthly  flows 
were  interpolated  in  a  similar  fashion. 
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4  Trends  and  Patterns  in  the 
Results 


Several  channel  processes  were  occurring  in  the  Gavins  Point  tailwater  that 
have  habitat  implications.  Degradation  of  the  channel  was  occurring  near  each 
of  the  dams  because  the  dam  disrupted  the  downstream  transport  of  sediment. 
The  degraded  sediment  was  transported  downstream  of  Omaha,  NE,  where  it 
was  deposited  in  the  aggrading  reaches  of  the  river.  These  observations  have 
been  made  on  numerous  occasions  and  are  not  unique  to  this  study;  however, 
when  viewed  from  the  perspective  of  an  RCHARC  analysis,  they  provide  a 
unique  insight  into  habitat  conditions  in  the  Missouri  River. 

The  habitat  changes  in  the  Missouri  River  are  best  described  by  comparing 
CSRS  (preproject)  transects  with  project  transects  at  similar  locations  (Fig¬ 
ures  12-15).  Several  significant  features  become  apparent.  First,  the  eleva¬ 
tions  in  the  project  channel  have  all  been  reduced  by  about  5  ft  fiom  the  CSRS 
channel.  Note  also  that  the  stage-discharge  is  in  general  compressed  in  the 
project  cross  sections  compared  to  the  CSRS  cross  sections.  The  compression 
is  greatest  for  the  two  transects  having  the  greatest  top  width  and  least  for  the 
two  transects  having  the  least  top  width.  Perhaps  the  lack  of  sediment  input  to 
this  reach  has  allowed  the  releases  from  the  dam  to  scour  out  a  more  hydrau¬ 
lically  efficient  channel.*  In  support  of  this  observation,  the  narrow  transects 
are  widening,  and  wider  transects  are  narrowing  when  tte  project  top  widths 
are  compared  to  the  preproject  system  as  though  the  river  is  “ditching.”  Ditch¬ 
ing  appears  to  be  occurring  in  the  wider  reaches  that  offer  the  greatest  habitat 
diversity  and  concomitantly  should  have  an  effect  on  the  fish  community 
(KaUemeyn  and  Novotny  1977). 

The  depths  associated  with  the  lowest  discharge  in  the  project  charmel  are 
greater  than  the  depths  associated  with  lowest  discharge  in  the  preproject  chan¬ 
nel  (Figures  12-15).  The  cause  of  the  change  in  depths  at  low  discharge 
between  the  preproject  and  postproject  channel  is  unknown.  However,  the 
deeper  initial  depths  are  compensated  by  the  reduced  increase  in  stage  as  dis¬ 
charge  increases  in  the  project  channel.  By  the  time  the  project  flows  reach 


^  Personal  Communication,  March  1992,  Brad  R.  Hall,  U.S.  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 
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Comparison  of  a  cross  section  (RM  778.9)  in  the  wide  category  collected  prior  to  impoundment  (left)  with  that  after  many  years 
of  regulation  (right).  Configuration  of  this  plot  is  the  same  as  in  Figure  12.  Note  that  the  project  elevations  have  been  reduced 
by  about  5  ft  by  degradation  and  that  the  stage-discharge  relationship  has  been  significantly  compressed 


50,000  cfs,  they  appear  to  have  much  the  same  depths  as  the  32,000-cfs  flows 
in  the  CSRS.  While  the  observations  in  the  previous  two  paragraphs  are  spec¬ 
ulative,  they  do  suggest  that  the  degradation  has  had  significant  habitat 
implications. 

Figures  16-19  present  plots  at  monthly  intervals  that  relate  the  similarity  as 
correlation  coefficients  of  the  historical  monthly  flows  to  range  of  incremental 
project  discharges.  Note  that  the  general  form  of  these  plots  follows  the  value 
function  format  required  by  MRD  to  link  the  “native  warmwater  riverine  fish 
community”  resource  category  to  the  other  resource  categories  as  part  of  the 
Missouri  River  Water  Control  Manual  Review  studies.  The  flows  used  for  the 
CSRS  are  listed  in  the  tabulation  in  the  section  “Hydrologic  Summaries”  under 
median  flows  for  Gavins  Point  Dam. 

Figures  16  and  17  present  top  width  unadjusted  correlation  coefficients  for 
depth  and  velocity,  respectively;  Figures  18  and  19  present  the  same  informa¬ 
tion  except  that  the  correlation  coefficients  have  been  adjusted  for  deviations 
from  historical  top  widths.  The  dotted  vertical  reference  line  represents  the 
median  flow  for  each  month.  The  solid  vertical  reference  line  on  each  of  the 
12  plots  represents  an  arbitrary  sequence  of  monthly  flows  selected  by  the 
authors  to  represent  an  imaginary  project  alternative. 

Several  trends  in  the  patterns  of  the  correlation  coefficients  are  apparent 
First,  the  correlation  coefficients  for  the  top  width  unadjusted  depth  correlation 
coefficients  have  an  apparent  bimodal  shape.  The  initial  minimum  in  the 
correlation  coefficients  for  each  month  represents  the  discrepancy  in  depth 
observed  at  the  lowest  simulated  flows.  The  first  peak  represents  the  point  at 
which  the  depth  distributions  of  the  CSRS  and  project  channels  converge. 

Depth  increases  less  quickly  with  discharge  in  the  project  channel  than  in  the 
CSRS.  Consequently,  the  depths  in  the  CSRS  “catch  up”  to  the  depths  in  the 
project  channel.  For  the  wide  channel  category  this  appears  to  happen  at  about 
20,000  cfs. 

The  correlation  coefficients  reach  a  minimum,  however,  when  the  dis¬ 
charges  in  the  CSRS  for  the  wide  and  divided  channel  categories  go  full  chan¬ 
nel  (full  channel  refers  to  point  at  which  the  stage  increases  enough  that  all  of 
the  channel  between  the  river  banks  is  inundated).  There  is  no  corresponding 
behavior  in  the  project  channel.  The  correlation  coefficients  then  again 
increase  as  the  flows  in  the  CSRS  go  higher  and  increase  their  depth  over  the 
shallow  parts  of  the  channel  and  become  deeper  like  the  project  channel.  The 
dip  in  the  columns  of  correlation  coefficients  corresponds  to  the  threshold 
when  the  historical  channel  acquires  a  high  percentage  of  shallow  depths  from 
going  full  bank  (Figures  12  and  15). 

Degradation  has  also  reduced  the  connectivity  of  the  main  channel  to  chutes 
and  backwaters  (Figure  15).  The  elevations  of  the  highest  discharges  are  inad¬ 
equate  to  flood  these  areas,  limiting  the  availability  of  shallower  water  at  high 
discharges.  Consequently,  side  channels  and  other  types  of  off-channel  habitat 
are  no  longer  a  habitat  feature  of  the  degrading  reaches  of  the  regulated 
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Figure  18.  Monthly  value  functions  (constructed  from  correlation  coefficients)  for  the  wide  channel  category  for 
depth  corrected  for  top  width.  The  conventions  of  Figure  16  are  followed  in  this  figure.  Note  that 
bimodal  peaks  can  be  observed  even  through  the  top  width  corrections 
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velocity  corrected  for  top  width.  The  conventions  of  Figure  16  are  followed  in  this  figure.  Note  that 
bimodal  peaks  can  be  observed  even  through  the  top  width  corrections 


Missouri  River  during  high  discharges.  Overbank  areas  are  not  available  for 
fishes,  and  organic  matter  produced  in  the  river  margins  no  longer  can  be 
inundated  by  the  flows.  Riverine  fishes  that  utilize  this  habitat  feature  could 
be  impacted  negatively. 

The  compression  of  the  stage-discharge  relationship  in  the  Gavins  Point 
tailwater  and  the  relatively  large  cross-sectional  area  available  at  lower  flows 
in  the  project  reach,  compared  to  the  preproject  reach,  also  distort  the  relation¬ 
ship  between  mean  cross-section  velocity  and  stage.  At  low  discharge,  the 
water  velocity  in  the  project  river  is  reduced  from  the  water  velocity  in  the 
CSRS.  However,  because  of  the  compression  of  the  stage-discharge  relation¬ 
ship  in  the  project  channel,  eventually  the  average  depth  of  the  project  reach 
will  be  the  same  as  that  of  the  CSRS  as  discharge  increases  in  both  channels. 
For  depth  in  the  wide  channel  category,  a  discharge  of  14,000  cfs  in  the  pro¬ 
ject  charmel  has  about  the  same  depth  distribution  as  a  discharge  of  25,000  cfs 
in  the  main  channel  (Figure  15). 

A  similar  pattern  of  disruption  would  occur  for  velocity  distributions;  that 
is,  the  discharge  at  which  the  velocity  distribution  between  the  project  and 
CSRS  channels  occurs  would  be  different  (and  also  different  from  the  dis¬ 
charges  at  which  the  depth  distributions  were  similar).  It  is  well-known  that 
depth  and  velocity  are  correlated  variables  in  rivers.  It  would  appear  that 
because  of  regulation  the  correlation  between  depth  and  velocity  is  different 
between  the  CSRS  and  the  project  channel.  While  these  conclusions  are  spec¬ 
ulative,  they  do  point  to  a  need  to  better  understand  how  regulation  may  have 
affected  the  physical  habitat  features  of  warmwater  fishes. 

The  change  in  the  depth/velocity  correlation  between  the  project  and  the 
CSRS  has  a  number  of  significant  ramifications  that  deserve  amplification. 

The  alterations  of  the  Missouri  River  from  its  preregulation  conditions  should 
be  studied  in  much  greater  detail  to  better  understand  the  effects  of  regulation 
and  as  a  prelude  for  any  restoration  efforts  on  the  Missouri  River.  The  cycles 
of  deposition  and  erosion  are  affected.  In  the  CSRS  channel,  the  low  flow  was 
more  of  an  erosional  feature  than  the  project  low-flow  charmel.  Conversely, 
because  the  present  high  flows  are  inadequate  to  reach  the  overbank  in  most 
years,  the  present  high  flows  that  are  relegated  to  the  main  charmel  are  more 
erosional  Aan  the  CSRS  high  flows  (less  conveyance  at  high  flows  in  the 
project  charmel)  and  less  likely  to  deposit  organic  matter  than  in  the  CSRS. 
Transport  and  storage  of  organic  matter  in  streams  is  well  known  to  affect  the 
composition  of  the  aquatic  community. 

Fish  habitat  is  affected  in  several  ways.  First,  because  the  correlation 
between  depth  and  velocity  has  been  altered  by  regulation,  formulation  of 
habitat  suitability  curves  for  detailed  population  studies  must  be  done  with 
considerable  caution.  The  present  habitat  provided  by  the  Gavins  Point  tail- 
water  is  considerably  different  than  provided  by  the  CSRS.  Suitability  curves 
developed  in  the  present  Missouri  River  may  be  misleading  if  the  effects  of 
charmel  alteration  are  not  considered.  Second,  the  availability  of  overbank 
areas  or  slower  velocity  areas  is  reduced  at  high  flows.  These  effects  are 
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probably  life  stage-specific  and  cannot  be  discussed  in  detail  in  this  report.  It 
is  almost  as  though  the  habitat  characteristics  associated  with  high  and  low 
flow  have  been  partially  reversed;  low  flows  in  the  project  channel  provide 
significant  cross-sectional  area  and  considerable  opportunity  for  storage  of 
orgamc  matter.  Missouri  River  in-channel  hydraulic  characteristics  are  now 
out  of  synchrony  with  the  hydraulic  patterns  in  the  CSRS. 

The  results  provided  by  the  RCHARC  analysis  are  summarized  in  Fig¬ 
ure  20.  Trends  and  patterns  in  the  correlations  can  be  easily  observed.  This 
can  be  important  in  exploring  alternatives  if  a  particular  alternative  provides 
good  habitat  (high  correlation  coefficients)  for  most  months  but  very  low  habi¬ 
tat  for  one  or  more  months.  In  Figure  20,  the  median  monthly  flow  is  pre¬ 
sented  as  the  holdout  flow.  Note  that  the  median  flow  does  not  provide  the 
maximum  habitat  value  (maximum  habitat  value  is  12)  because  of  channel 
changes  and  distortions  of  the  stage-discharge  relationship. 
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Figure  20.  Summary  results  for  an  application  of  the  RCHARC  to  the  Gavins  Point  tailwater. 
The  upper  left  plot  provides  an  estimate  of  the  median,  25  percent  exceedance 
flow,  and  75  percent  exceedance  flow  for  the  CSRS.  An  arbitrary  annual  flow 
sequence  of  monthly  flows  represents  the  project  flows.  The  flows  are  presented 
on  a  log  scale  because  the  natural  logarithm  of  discharge  is  linearly  related  to 
depth.  The  log  scale  gives  an  approximation  of  how  depth  could  change  as  flows 
change.  The  upper  right  plot,  the  deviation  plot,  describes  the  degree  to  which  the 
flows  differ  from  the  median.  The  lower  two  plots  represent  an  annual  time  series 
of  monthly  correlation  coefficients  associated  with  the  holdout  (in  this  case,  the 
median  flow)  and  the  operational  flow.  The  sum  of  the  holdout  correlation  coeffi¬ 
cients  represents  the  effects  of  the  median  CSRS  hydrograph;  the  sum  of  the  oper¬ 
ational  correlation  coefficients  represents  the  impacts  of  the  project  flow 
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5  Sources  of  Error  and 
Uncertainty 


As  with  any  new  methodology,  RCHARC  has  produced  the  following  sug¬ 
gestions  for  improving  the  approach,  sources  of  error,  and  important  caveats 
and  assumptions: 

a.  The  discharges  at  which  the  CSRS  cross  sections  were  taken  are 
unknown.  The  degree  to  which  the  cross  sections  may  have  changed  at 
other  discharges  is  unknown.  Because  the  two  calibration  cross  sec¬ 
tions  collected  at  each  site  for  the  project  data  differed  little,  this  is  not 
considered  a  major  problem. 

b.  It  is  not  possible  to  detennine  the  error  in  the  velocity  predictions. 
Although  the  cross  sections  are  reasonably  similar  between  the  CSRS 
data  set  and  the  project  channel,  it  was  not  learned  that  the  slope  of  the 
stage-discharge  relationship  varies  between  CSRS  and  project  channels 
until  the  study  was  almost  complete.  The  differences  in  relative  stage 
between  the  two  data  sets  probably  introduced  a  bias  into  the  quadratic 
equation  used  for  the  velocity  correction.  Also,  velocity  is  a  derived 
variable  (based  on  a  calibration  data  set  and  not  directly  predicted  like 
depth)  in  the  project  chaimel  and  a  twice-derived  variable  in  the  CSRS. 
In  the  CSRS,  velocities  are  first  synthesized  based  on  hydraulic  radius 
and  then  adjusted  with  a  correction  factor. 

c.  As  with  any  method  of  smoothing,  some  loss  of  information  occurs 
when  curves  are  employed  instead  of  original  data.  However,  given  the 
highly  significant  fits  of  the  curves  to  the  data,  information  loss  associ¬ 
ated  with  curve  fitting  is  not  considered  to  be  a  problem,  and  in  fact, 
curve-fitted  data  probably  more  accurately  capmre  trends  and  patterns 
of  depth  and  velocity  in  the  transects. 

d.  Considerable  problems  were  experienced  in  obtaining  a  consistent 
datum  for  the  field  data.  Although  good  hydrologic  practice  was  fol¬ 
lowed,  it  is  possible  that  stage  errors  remain  in  the  analysis. 
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e.  The  method  used  to  correct  the  velocity  predictions  for  the  CSRS  chan¬ 
nel  will  result  in  spurious  self-correlation  when  the  CSRS  velocity 
distributions  are  compared  to  the  project  channels  because  a  proportion 
of  the  distribution  pattern  in  the  project  channel  (that  part  represented 
by  the  quadratic  equation)  also  occurs  in  the  CSRS.  Thus,  die  correla¬ 
tion  coefficients  at  equivalent  discharges  will  tend  to  be  inflated.  While 
the  velocity  correlations  are  useful  and  can  provide  insight  into  habitat 
conditions  in  the  Missouri  River,  they  should  not  receive  the  same 
credibility  as  the  depth  correlation  coefficients  unless  sensitivity  analy¬ 
ses  are  performed  with  both  depth  and  velocity  to  deteraiine  if  the 
pattern  in  the  results  is  consistent  between  the  two  variables. 

/.  Multiple  sources  of  existing  and  newly  collected  hydraulic  and  hydro- 
logic  data  were  employed  in  this  study.  There  was  no  opportunity  to 
independently  verify  some  of  the  data  in  this  study,  particularly  for  the 
CSRS  charmel. 

g.  Impact  will  be  expressed  as  a  single  number  that  integrates  impact  over 
a  93-year  hydrologic  record.  However,  there  are  many  different  pat¬ 
terns  that  will  generate  the  same  final  estimate  of  total  impact.  The 
sequence  of  correlation  coefficients  should  be  evaluated  to  separate  the 
alternatives  that  generate  many  years  of  average  habitat  from  those 
alternatives  that  generate  a  mixture  of  extremely  high  and  extremely 
low  habitat  values. 
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6  Discussion 


The  similarity  results  provided  by  the  RCHARC  analysis  are  shown  in 
Appendix  B.  In  the  same  manner  as  in  the  univariate  analysis  using  correla¬ 
tions,  the  Canberra  coefficients  can  be  used  to  assess  the  various  alternatives  to 
distinguish  the  conditions  most  similar  to  historical  hydraulic  patterns  based  on 
median  monthly  flows.  Seasonally  dependent  channel  habitat  conditions  most 
like  those  of  the  preproject  era  can  be  approached  on  a  species  or  life  history- 
specific  basis. 

Tables  of  Canberra  coefficients  such  as  those  in  Appendix  B  give  resource 
managers  an  opportumty  to  choose  alternatives  based  on  similarity  matrices  for 
a  given  water  year  (in  this  case  a  median  year).  However,  from  these  numbers 
it  is  impossible  to  identify  actual  changes  in  depth-velocity  profiles  from  his¬ 
torical  to  existing  conditions.  To  allow  for  visual  observation  of  hydraulic 
habitat  alterations,  contour  plots  (Figures  21  and  22)  were  generated  that  illus¬ 
trate  preproject  and  project  depth-velocity  relationships  relative  to  total 
proportions. 

Contour  plots  of  two  types  were  produced:  (a)  plotting  depth  versus  veloc¬ 
ity  with  the  function  variable  percent  of  total  frequency  (per  discharge)  used  as 
the  contouring  parameter  (for  historical  and  existing  conditions,  respectively); 
and  (b)  plotting  depth  versus  velocity  with  the  function  variable  percent  fre¬ 
quency  difference  (the  difference  between  historical  and  project  frequency 
distributions)  used  as  the  contouring  parameter.  The  latter  plots  present  easily 
discernible  information  regarding  loss  and  gain  in  hydraulic  habitat  from  pre¬ 
project  to  project  conditions.  To  compliment  the  habitat  quality  infonnation, 
quantity  was  included  in  the  form  of  top  width  within  each  plot. 

Figure  21  depicts  narrow  channel  plots  for  flows  in  a  median  winter  month 
(10,000  cfs  and  28,000  cfs)  for  preproject  and  project  conditions,  respectively. 
The  contour  plot  of  historic  conditions  reveals  a  diverse  set  of  available  habitat 
types  including  deep-slow  and  deep-fast,  with  higher  proportions  of  habitat 
with  moderate  depth  and  slow-moving  water.  Examination  of  existing  depth- 
velocity  conditions  indicates  a  more  uniform  and  less  diverse  range  of  avail¬ 
able  habitats,  yielding  a  direct  relationship  between  the  two  variables  (velocity 
increases  with  depth — a  relationship  not  apparent  in  the  fonner  plot).  Also 
note  the  complete  absence  of  deep-slow  water  and  the  addition  in  high 
proportion  of  deeper-faster  water  areas.  The  bottom  plot  in  Figure  21 
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Figure  21 .  Bivariate  contour  plots  for  narrow  channel.  Top  and  middle  plots  relate  depth  and 
velocity  (as  a  bivariable)  to  proportion  of  total  frequency  each  depth-velocity  incre¬ 
ment  represents.  Upper  plot  presents  historical  conditions;  middle  plot  presents 
existing  conditions.  Discharges  modeled  represent  median  winter  monthly  flows 
(10,000  cfs  for  preproject;  28,000  cfs  for  project).  Note  the  more  diverse  habitat 
conditions  available  in  historical  plot  and  more  uniform,  direct  depth-velocity  rela¬ 
tionship  evident  in  project  plot.  Quantity  of  available  habitat  is  included  in  these 
two  plots  in  the  form  of  top  width.  Bottom  plot  shows  the  differences  in  percent 
frequency  between  historical  and  existing  depth-velocity  conditions.  Contour  lines 
depict  gains  in  habitat  and  lighter  broken  lines  depict  habitat  loss.  Top  diff 
(preproject-project  top  width)  signifies  change  in  habitat  quantity.  Contour  levels 
are  shown  in  the  legend 
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Figure  22.  Preproject  contour  plots  and  difference  in  frequency  plots  for  wide  and  transitional 
channels.  Median  monthly  winter  flows  of  28,000  cfs  were  used  in  modeling  project 
conditions  for  both  channels.  Note  the  almost  complete  loss  of  originally  available 
habitat  and  general  gain  in  habitats  of  greater  depths  in  both  frequency  difference 
plots 
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illustrates  the  loss  and  gain  of  habitat  type  in  a  narrow  channel  for  these 
particular  flows. 

Changes  in  depth-velocity  relationships  from  preproject  to  project  condi¬ 
tions  for  wide  and  transitional  channels  are  shown  in  Figure  22.  The  project 
contour  plots  are  omitted  here,  assuming  the  pertinent  information  can  be 
gleaned  from  the  difference  in  frequency  plots.  As  in  Figure  21,  these  contour 
plots  represent  historical  flows  of  10,000  cfs  and  existing  flows  of  28,000  cfs, 
and  the  top  difference  (diff)  reflects  the  project  top  width  subtracted  from  the 
preproject  top  width.  Preproject  conditions  for  the  wide  channel  indicate  habi¬ 
tats  of  mostly  slow-moving  shallow  water  with  some  swifter  water  at  shallow 
and  moderate  depths.  The  difference  in  frequency  plot  for  this  channel  reveals 
a  loss  of  most  of  the  originally  available  habitat  and  the  addition  of  two  pre¬ 
viously  unavailable  habitat  types;  slow,  shallow  water  and  swift  water  at 
greater  depths.  For  the  transitional  channel,  preproject  conditions  show  a 
dominance  of  shallow,  mostly  slow  water  for  the  median  monthly  winter  flow. 
The  difference  in  frequency  plot  for  this  channel  indicates  that  with  increased 
flows  and  a  more  narrow  top  width,  habitats  of  swifter  water  with  considerable 
depth  gradients  become  available.  The  slow,  shallow  water  apparent  in  the 
historical  conditions  is  absent  in  the  existing  conditions. 
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7  Conclusions  and 
Recommendations 


The  study  generated  the  following  conclusions  and  recommendations; 

a.  The  correlation  coefficients  generated  for  the  Gavins  Point  tailwater 
could  be  interpreted  and  explained;  in  fact,  correlation  may  be  too 
sensitive  a  statistic  for  all  applications  of  the  RCHARC.  The  RCHARC 
provided  insight  to  trends  and  patterns  of  physical  habitat  in  this  reach 
that  could  not  have  been  determined  in  any  other  way.  Based  on  the 
Gavins  Point  analysis,  the  RCHARC  was  able  to  summarize  the  effects 
of  different  operational  alternatives  on  habitat  for  warmwater  fishes, 
assuming  that  the  CSRS  represented  the  habitat  ideal  for  this 
community. 

b.  Given  the  multiple  uncertainties  in  the  velocity  data  sets,  use  of  the 
depth  correlation  coefficients  as  the  basis  of  the  impact  assessment  is 
suggested.  The  velocity  correlation  coefficients  can  be  used  if  sensi¬ 
tivity  analysis  suggests  that  the  two  variables  are  providing  much  the 
same  information. 

c.  Velocity  distributions  should  be  considered  and  discussed  as  informa¬ 
tion  supplemental  to  the  depth  information. 

d.  The  Gavins  Point  tailwater  of  the  Missouri  River  is  often  presented  as 
the  last  remaining  “natural”  reach  in  the  system.  This  analysis  has 
described  a  number  of  habitat  alterations  resulting  from  regulation  in 
this  reach.  A  truly  “natural”  segment  of  the  lower  Missouri  River  may 
no  longer  exist. 

e.  The  RCHARC  analysis  both  quantified  and  explained  alterations  in  the 
Gavins  Point  reach  of  the  river.  An  expansion  of  the  RCHARC  analy¬ 
sis  would  be  valuable  for  restoration  planning  on  the  main  stem 
Missouri  River  or  any  other  large  river  considered  for  restoration. 
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/.  The  RCHARC  analysis  was  performed  at  an  ecosystem  level  to 

describe  major  changes  in  the  hydraulics  and  hydrology  of  the  Gavins 
Point  reach.  The  patterns  and  trends  identified  in  this  analysis  should 
focus  on  future  fishery  and  ecological  investigations. 

g.  The  RCHARC  analysis  considered  only  physical  habitat  changes 
between  the  CSRS  and  the  project  channel.  Qearly,  habitat  changes 
are  a  major  feature  of  the  impacts  of  different  alternatives;  however, 
because  hydraulics  affects  everything  else,  other  variables  such  as  water 
temperature  changes,  water  turbidity,  and  nutrient  cycles  should  also  be 
evaluated. 

h.  Statistics  other  than  the  median  could  be  used  to  summarize  the  long¬ 
term  impact  of  each  operational  alternative.  Possible  options  include 
minimum  habitat  analysis,  increased  weighting  of  key  months,  and  the 
use  of  special  rules  (two  sequential  poor  habitat  years  are  emphasized). 

i.  Like  many  analysis  methodologies,  the  RCHARC  can  generate  numbers 
and  the  numbers  can  be  ranked,  summed,  and  processed  to  support 
conclusions.  However,  it  is  critical  that  a  RCHARC  analysis  be  done  at 
a  level  of  rigor  to  understand  differences  in  fluvial  processes  between 
the  CSRS  and  the  project  alternatives  that  have  ecological  significance. 

J.  Habitat  for  threatened  and  endangered  species  needs  to  be  considered  in 
light  of  RCHARC  analysis  first  and  then  in  context  of  standard  habitat- 
based  approaches.  Development  of  habitat  suitabilities  for  threatened 
and  endangered  species  must  acknowledge  the  major  shifts  in  the  depth 
and  velocity  patterns  that  have  occurred  in  the  Missouri  River,  even  in 
the  reach  that  is  normally  considered  “natural.”  Certain  critical  habitat 
features  may  be  much  less  abundant  in  the  project  river.  It  is,  there¬ 
fore,  difficult  to  determine  utilization  that  allows  complete  development 
of  habitat  suitability  curves. 

k.  The  interpretation  of  the  RCHARC  analysis  should  also  be  performed 
using  a  bivariate  approach  for  the  Fort  Randall,  Garrison,  and  Fort  Peck 
tailwaters  to  at  least  the  rigor  of  the  Gavins  Point  analysis  to  describe 
and  understand  the  effects  of  river  regulation  on  downstream  habitat. 
The  effects  of  river  regulation  on  the  warmwater  fish  community  must 
be  fully  understood,  particularly  for  threatened  and  endangered  species 
whose  decline  has  been  associated  with  regulation. 

/.  Bivariate  contour  plots  based  on  proportional  frequencies  present  easily 
discernible  information  regarding  presence  and  absence  of  depth- 
velocity  profiles  through  space  (chaimel  category)  and  time  (monthly, 
yearly,  or  any  other  applicable  time  scale).  These  plots  allow  resource 
managers  rapid  visual  assessment  of  available  depth-velocity  contours 
given  particular  objectives  (i.e.  recreational  uses,  fishery  management) 
and  flow  alternatives. 
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The  bivariate  contour  plots  can  also  aid  biologists  and  ecologists  con¬ 
cerned  with  threatened  and  endangered  species  management,  impact 
assessment  and  restoration,  and  community  level  management.  Habitat 
preferences  for  endangered  species,  fish  communities,  or  guilds  (i.e. 
riffle  dwellers)  can  be  superimposed  onto  the  contour  plots  to  deteraiine 
if  habitat  or  life  history  requirements  are  achieved.  The  biologist  can 
quickly  peruse  all  possible  flow  alternatives  and  decide  which  flow 
condition  produces  the  best  patterns  of  habitat. 
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Table  1 

Example  of  IFG4  Output  Processed  by  the  LSTVDX  Program 
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0.00 

90.00 

B 

2043.000 

2023.000 

0.00 

0.00 

90.00 

B 
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2.65 

1.01 

90.00 

B 

2087.000 

2067.000 

9.15 

2.31 

90.00 

B 

2107.000 

2087.000 

7.65 

2.05 

90.00 

B 

2127.000 

2107.00 

10.35 

2.51 
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B 
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2127.00 

14.85 

3.19 
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B 
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3.64 
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B 
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17.15 

3.52 

90.00 

B 
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17.65 
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3.61 
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4.85 

1.51 
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2557.000 

2517.000 

4.35 

1.40 
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2607.000 

2557.000 
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0.98 
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2657.000 
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0.26 
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1.01 
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1.13 
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1.01 
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1.22 
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5.35 

1.61 
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1.96 
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B 
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8.55 

2.21 
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5.05 

1.55 
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3207.000 

0.00 

0.00 
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3240.000 
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0.00 
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90.00 

B 

3267.000 
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0.00 
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Note: 

Row  A  contains  the  river  mile  discharge,  and  water  surface  elevation.  Row  B  contains  the  left 

and  right  bank  adjacent  cell  boundaries,  cell  depth,  cell  velocity,  and  cell  cover. 

1 

Table  2 

Summary  Statistics  for  Logistics  Equation  Fitted  to  the 

Logistics  Equation  for  the  Cumulative  Depth  Frequency 
Distribution 

-  -  —  r.MAMnATF  niv/incn  nicrUADriC  cnnn 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

265388.36125  88462.78708 

Residual 

33 

196.13955  5.94362 

Uncorrected  Total 

36 

265584.50080 

(Corrected  Total) 

35 

16997.21240 

CHANCATF 

ni\/incn  ruc^UAnr^c:  orsr\r\ 

Source 

DF 

ovA^u - 

Sum  of  Squares  Mean  Square 

Regression 

3 

276241.42995  92080.47665 

Residual 

35 

143.39165  4.09690 

Uncorrected  Total 

38 

276384.82160 

(Corrected  Total) 

37 

23361.04337 

CHANCATF 

nix/incn  rMor'UAnoc:  li 

Source 

DF 

■  1  tvwu - 

Sum  of  Squares  Mean  Square 

Regression 

3 

285582.97976  95194.32659 

Residual 

36 

194.31627  5.39767 

Uncorrected  Total 

39 

285777.29603 

(Corrected  Total) 

38 

24209.77055 

--  -  __  CHANCATE 

~rii\/incn  nioruADric  4/iAAn  1 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

284169.10381  94723.03460 

Residual 

36 

204.11066  5.66974 

Uncorrected  Total 

39 

284373.21447 

(Corrected  Total) 

38 

24725.08661 

CHANCATE 

-r»i\/jnpn  nioruADnc  -icaaa 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

292358.03272  97452.67757 

Residual 

37 

165.26034  4.46650 

Uncorrected  Total 

40 

292523.29306 

(Corrected  Total) 

39 

26056.44776 

CHANCATE' 

-nivinpn  nicruADnc  oaaaa 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

236273.03927  78757.67976 

Residual 

32 

353.54431  11.04826 

Uncorrected  Total 

35 

236626.58358 

(Corrected  Total) 

34 

25577.73558 
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1  Note: 

Note  that  the  F-statistic  (obtained  by  dividing  the  Mean  Square  Error  of  the  Regression  by  the 

Mean  Square  of  the  Residuals)  varies  between  800  and  10,000.  This  range  of  F-values  for 

linear  regression  would  be  highly  significant. 
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CHANCATE=WIDE  DISCHARGE=24000 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

250844.80073 

83614.93358 

Residual 

35 

186.37862 

5.32510 

Uncorrected  Total 

38 

251031.17935 

(Corrected  Total) 

37 

25351.37738 

CHANCATE=WIDE  DISCHARGE=28000 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

264462.24240 

88154.08080 

Residual 

36 

332.70368 

9  .24177 

Uncorrected  Total 

39 

264794.94608 

(Corrected  Total) 

38 

23525.46113 

CHANCATE=WIDE  DISCHARGE=32000 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

272801.66684 

90933.88895 

Residual 

37 

462.43630 

12.49828 

Uncorrected  Total 

40 

273264.10314 

(Corrected  Total) 

39 

25365.87897 

CHANCATE=WIDE  DISCHARGE=36000 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

298889.60924 

99629.86974 

Residual 

39 

502.76844 

12.89150 

Uncorrected  Total 

42 

299392.37768 

(Corrected  Total) 

41 

24277.681 15 

-  CHANCATE=WIDE  DISCHARGE=40000 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

301489.90375 

100496.63458 

Residual 

40 

472.75769 

11.81894 

Uncorrected  Total 

43 

301962.66144 

(Corrected  Total) 

42 

28722.56358 

-  CHANCATE=WIDE  DISCHARGE=46000 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

314097.73009 

104699.24366 

Residual 

42 

346.39820 

8.24758 

Uncorrected  Total 

45 

314444.12829 

(Corrected  Total) 

44 

32042.78763 

-  CHANCATE=WIDE  DISCHARGE=50000 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

309407.11366 

103135.70455 

Residual 

42 

341.75651 

8.13706 

Uncorrected  Total 

45 

309748.87017 

(Corrected  Total) 

44 

34715.36237 
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-  CHANCATE=WIDE  DISCHARGE=6000  - 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

164914.15726 

54971.38575 

Residual 

21 

201.62280 

9.60109 

Uncorrected  Total 

24 

165115.78006 

(Corrected  Total) 

23 

17778.79290 

OU  AKir' ATC. 

wir 

niC5nHARttF=P 

JOOO  -  -  -  - 

- OnMlMOM  1  C' 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

157294.08089 

52431.36030 

Residual 

21 

273.01540 

13.00073 

Uncorrected  Total 

24 

157567.09629 

(Corrected  Total) 

23 

20932.23726 

r'U  AMP  ATC 

_\A/ir 

niQr.uADnc-' 

1 0000  - 

-■  - —  OnMINvM  1  c 

Source 

DF  Sum  of  Squares 

Mean  Square 

Regression 

3 

158604.68939 

52868.22980 

Residual 

22 

238.49590 

10.84072 

Uncorrected  Total 

25 

158843.18259 

(Corrected  Total) 

24 

24254.28964 

CHANCATE=WIDE  DISCHARGE=12000 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

163919.29259 

54639.76420 

Residual 

23 

180.74781 

7.85860 

Uncorrected  Total 

26 

164100.04040 

(Corrected  Total) 

25 

25496.61848 

PU AMP ATC 

-wir 

\c  ntcpuADriP-iAnnn - 

- LTrlMINs/M  1  C 

Source 

—  VVIL 

DF 

1- 

Sum  of  Squares 

Mean  Square 

Regression 

3 

158572.55655 

52857.51885 

Residual 

23 

172.33217 

7.49270 

Uncorrected  Total 

26 

158744.88872 

(Corrected  Total) 

25 

26458.21832 

pi_iAMpATc_\A/tnc  mcpuADrsP-ifinnn - - 

- 1  ini  i  > 

Source  DF  Sum  of  Squares 

Mean  Square 

Regression 

3 

165361.09117 

55120.36372 

Residual 

24 

128.65409 

5.36059 

Uncorrected  Total 

27 

165489.74526 

(Corrected  Total) 

26 

27189.91800 

CHANCATE=WIDE  D1SCHARGE=20000 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

235232.07482 

78410.69161 

Residual 

34 

180.30950 

5.30322 

Uncorrected  Total 

37 

235412.38432 

(Corrected  Total) 

36 

30957.34345 
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/^UAKinATC 

TDAMQiT  nicruAnnF  OAnnn  —  —  .  —  —  - 

---■  - - - — 1  ^ 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

148884.19251  49628.06417 

Residual 

29 

188.95754  6.51578 

Uncorrected  Total 

32 

265584.50080 

(Corrected  Total) 

31 

16997.21240 

/^UAMr'ATP 

-TDAMQIT  niQP.UADriD  OPHnO  _  _  -  —  - 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

146261.03154  48753.67718 

Residual 

29 

193.96908  6.68859 

Uncorrected  Total 

32 

146455.00061 

(Corrected  Total) 

31 

36668.21689 

r'UIAKir'ATP 

TR/ 

niSHHARfiF  32000  - 

— -  - L/riMINwM  1  C 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

1 53487.96026  51 1 62.65342 

Residual 

30 

164.71370  5.49046 

Uncorrected  Total 

33 

153652.67397 

(Corrected  Total) 

32 

34855.71577 

rWANPATF 

TDAMQIT  niQr.UADriD  QRono 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

159722.41820  53240.80607 

Residual 

32 

166.16008  5.19250 

Uncorrected  Total 

35 

159888.57828 

(Corrected  Total) 

34 

37650.67619 

r'UAM^ATC 

TDAMQIT  niQr.uAPriF  Annnrs . 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

164031.54347  54677.18116 

Residual 

33 

205.89148  6.23914 

Uncorrected  Total 

36 

164237.43495 

(Corrected  Total) 

35 

36855.30653 

OUAKir'ATC 

TDAMQIT  niQruAoriF  iipmn  .  - 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

173455.48732  57818.49577 

Residual 

35 

172.68530  4.93387 

Uncorrected  Total 

38 

173628.17261 

(Corrected  Total) 

37 

37634.01262 

PHANrATF 

TDAMQIT  niQruiADriF  - 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

176623.94008  58874.64669 

Residual 

36 

198.76118  5.52114 

Uncorrected  Total 

39 

176822.70126 

(Corrected  Total) 

38 

38524.19515 

. 
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Table  2  (Continued) 


-  CHANCATE 

Source 

E=DI\ 

DF 

nOED  DISCHARGE: 

Sum  of  Squares 

=24000  - 

Mean  Square 

Regression 

3 

247519.30844 

82506.43615 

Residual 

34 

244.80791 

7.20023 

Uncorrected  Total 

37 

247764.11636 

(Corrected  Total) 

36 

30334.70309 

CHANCATE 

-oonnn 

Source 

DF 

•  ir-ii 

Sum  of  Squares 

~tOV/Uv  —  — 

Mean  Square 

Regression 

3 

245734.08387 

81911.36129 

Residual 

34 

218.67071 

6.43149 

Uncorrected  Total 

37 

245952.75458 

(Corrected  Total) 

36 

30077.09131 

CHANCATE 

:^ni\/irtcn  nfoouADr'C-ooAAA 

v/i  1  c 

- 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

256871.73201 

85623.91067 

Residual 

35 

175.91347 

5.02610 

Uncorrected  Total 

38 

257047.64548 

(Corrected  Total) 

37 

29273.43862 

CHANCATE 

-rti\/inpn  rtiQr.uADric_ocnAn 

Source 

DF 

Sum  of  Squares 

-wwwv/ - - - 

Mean  Square 

Regression 

3 

260060.00144 

86686.66715 

Residual 

36 

170.70553 

4.74182 

Uncorrected  Total 

39 

260230.70698 

(Corrected  Total) 

38 

32843.17545 

CHANCATE 

~rti\/fncrt  rticruADric_iinnnn 

Source 

DF 

Sum  of  Squares 

- 

Mean  Square 

Regression 

3 

263398.33931 

87799.44644 

Residual 

37 

154.50805 

4.17589 

Uncorrected  Total 

40 

263552.84736 

(Corrected  Total) 

39 

35582.29413 

CHANCATE 

—  •-✓I  V  1 

HUUMW  - 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

256348.85239 

85449.61746 

Residual 

37 

141.02930 

3.81160 

Uncorrected  Total 

40 

256489.88169 

(Corrected  Total) 

39 

37989.78081 

CHANCATE=DIVIDED  DISCHARGE=50000 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

259684.67020 

86561.55673 

Residual 

38 

159.33291 

4.19297 

Uncorrected  Total 

41 

259884.00311 

(Corrected  Total) 

40 

41816.06876 

Table  2  (Continued) 

rHAMP.ATF 

MADDOXA/  nicouADr^c^cnnn 

- -  “  -  — — 1  c 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

180344.36916 

60114.78972 

Residual 

31 

2398.68840 

77.37705 

Uncorrected  Total 

34 

182743.05756 

(Corrected  Total) 

33 

47985.19447 

--  _  _  PHANPATF 

MADorxvA/  npcruADr'c 

firmn 

— OvA/U 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

180029.60872 

60009.86957 

Residual 

32 

2713.44884 

84.79528 

Uncorrected  Total 

35 

182743.05756 

(Corrected  Total) 

34 

51835.41913 

_  PHAMPATF 

MADDnvA/  nioruADnc 

innnn 

Source  D 

F 

Sum  of  Squares 

Mean  Square 

Regression 

3 

176364.491 18  58788.16373 

Residual 

32 

2819.36468 

88.10515 

Uncorrected  Total 

35 

179183.85586 

(Corrected  Total) 

34 

52433.37459 

_  PHAMPATF 

KiADDrxXA/  nicnuADr^c 

lonnn 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

176940.38899 

58980.12966 

Residual 

33 

3263.81681 

98.90354 

Uncorrected  Total 

36 

180204.20580 

(Corrected  Total) 

35 

54348.19673 

_  PHAISIPATF 

-MAODnXA/  ntOnUADriC 

1  iionn 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

176638.63821 

58879.54607 

Residual 

34 

3565.56760 

104.86964 

Uncorrected  Total 

37 

180204.20580 

(Corrected  Total) 

36 

57749.71049 

-  _  PHANPATF 

MADDOVA/  nicruADr^c 

iRnrm  - 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

173915.71599 

57971 .90533 

Residua! 

34 

3381.12676 

99.44490 

Uncon'ected  Total 

37 

177296.84275 

(Corrected  Total) 

36 

57900.14660 

_ CHAWPATF- 

-MADDrxVA/  niQPUADr^r 

-ooooo  _ 

wn^i  1 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

134544.22265 

44848.07422 

Residual 

27 

322.30776 

11.93732 

Uncorrected  Total 

30 

134866.53041 

(Corrected  Total) 

29 

31305.30231 
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Table  2  (Continued) 

PHAMPATP 

NAP 

5RnW  ni9PHAR(5F 

24000  _  _  -  --  — 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

139257.65773 

46419.21924 

Residua! 

28 

238.66998 

8.52393 

Uncorrected  Total 

31 

139496.32771 

(Corrected  Total) 

30 

31392.26413 

PHANPATF 

-MADDOVA/  niCpUAPPC: 

28000  _  ~  - 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

149427.90958 

49809.30319 

Residual 

30 

161.62847 

5.38762 

Uncorrected  Total 

33 

149589.53804 

(Corrected  Total) 

32 

34261.94473 

--  PHAMPATF 

-MADDnvA/  nior.uADr^c: 

32000 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

145123,08041 

48374.36014 

Residual 

31 

157.55730 

5.08249 

Uncorrected  Total 

34 

145280.63771 

(Corrected  Total) 

33 

36702.77300 

ruAMrATTc  MADDriVA/  nicruADrsr 

qcnno 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

151618.34925 

50539.44975 

Residual 

33 

150.33897 

4.55573 

Uncorrected  Total 

36 

151768.68823 

(Corrected  Total) 

35 

41937,69929 

PHAMPATF 

MADDOVA/  nioruAD/^c  iinnno 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

163275.05655 

54425.01885 

Residual 

36 

142.43654 

3.95657 

Uncorrected  Total 

39 

163417.49309 

(Corrected  Total) 

38 

46660.14727 

PHAMPATF 

MADDnVA/  nicruAD/^c 

“46000 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

182259.00020 

60753.00007 

Residual 

37 

299.82518 

8.10338 

Uncorrected  Total 

40 

182558.82538 

(Corrected  Total) 

39 

37804.43083 

_  _  _  _  PHAMPATF 

-KJADOnW  niQPUAPPF 

50000  -  - 

Source 

DF 

Sum  of  Squares 

Mean  Square 

Regression 

3 

183610.51235 

61203.50412 

Residual 

38 

371.07824 

9.76522 

Uncorrected  Total 

41 

183981.59059 

(Corrected  Total) 

40 

42019.87349 
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Table  2  (Concluded) 

rMANrATF 

TDAMoiT  ntoruADric:  cnm  _  _  __  _ 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

5301 3.350870  1 7671 . 1 1 6957 

Residual 

11 

50.260730  4.569157 

Uncorrected  Total 

14 

53063.611601 

(Corrected  Total) 

13 

14209.541080 

CHANCATE 

-TDAWCIT  niQCUADCF-Qnnn  ....  -  _  i 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

56844. 1 56635  1 8948 .0522 1 2 

Residual 

12 

52.122412  4.343534 

Uncorrected  Total 

15 

56896.279047 

(Corrected  Total) 

14 

15245.577630 

CHANCATE 

-TDAMQiT  niQPUAPCF-innnn _ 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

601 79.623312  20059.874437 

Residual 

13 

63.876907  4.913608 

Uncorrected  Total 

16 

60243.500220 

(Corrected  Total) 

15 

15652.239235 

CHANPATE 

TDAMQIT  niQruADriC  lonnn 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

6481 6.485981  21 605.495327 

Residual 

14 

99.050356  7.075025 

Uncorrected  Total 

17 

64915.536336 

(Corrected  Total) 

16 

15330.741086 

CHANCATE 

-TDAMQIT  niQPUADriD-iiinnn  _  _ 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

69063.21 6803  23021 .072268 

Residual 

15 

183.892873  12.259525 

Uncorrected  Total 

18 

69247.109676 

(Corrected  Total) 

17 

17359.067286 

_ _  CHAMPATF 

-TDAMQIT  niQnuADnc  icnnn  1 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

69874.573443  23291 .524481 

Residual 

15 

376.517068  25.101138 

Uncorrected  Total 

18 

70251.090511 

(Corrected  Total) 

17 

14283.250243 

_  CHANCATE 

-TDAMQIT  niQPUADCF-onnnn  _  _ 

_ 

Source 

DF 

Sum  of  Squares  Mean  Square 

Regression 

3 

147681 .27493  49227.09164 

Residual 

28 

157.82264  5.63652 

Uncorrected  Total 

31 

147839.09757 

(Corrected  Total) 

30 

33046.05280 
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Table  3 


Appendix  A 

Top  Widths  and  Correiation 
Coefficients  of  the  Gavins 
Point  Taii  water 


Tables  of  coefficients  are  presented  in  two  major  groups.  The  first  group 
contains  the  correlation  coefficients  for  depth,  and  the  second  contains  the 
correlations  for  velocity.  Within  each  group,  the  tables  occur  in  three  major 
subsets,  one  subset  for  each  of  the  types  of  water  year  (median  flow,  high 
flow,  and  low  flow).  Within  each  subset  the  coefficients  are  separated  by 
channel  category  and  then  by  top  width  adjustment.  Within  each  pair  of 
tables,  the  first  table  contains  coefficients  not  adjusted  for  top  width  and  the 
second  member  of  the  pair  contains  coefficients  adjusted  for  top  width. 


GAVINS  CHANNEL  TOPWIDTHS  BY  DISCHARGE  AND  CHANNEL  CATEGORY 


OBS 

Q 

DIVOTOP 

NAROPTOP 

TRAOPTOP 

WIDOPTOP 

1 

6000 

1237 

1185 

909 

1939 

2 

8000 

1237 

1185 

981 

1974 

3 

10000 

1391 

1185 

1055 

2082 

4 

12000 

1678 

1204 

1148 

2107 

5 

14000 

1678 

1204 

1192 

2207 

6 

16000 

1768 

1204 

1377 

2275 

7 

20000 

2056 

1228 

1425 

2001 

8 

24000 

2198 

1322 

1457 

2397 

9 

28000 

2294 

1385 

1513 

2758 

10 

32000 

2528 

1428 

1578 

2956 

11 

36000 

2638 

1441 

1672 

3455 

12 

40000 

2756 

1496 

1822 

3549 
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GAVINS  UNADJUSTED  CORRELATION  COEFFICIENTS  OF  HISTORICAL  AND 
OPERATIONAL  DEPTHS  FOR  DIVIDED  CHANNEL  CATEGORY  -  AVERAGE  YEAR 
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Appendix  B 

Gavins  Point  Canberra 

Coefficients 


Tables  of  coefficients  of  historical/operational  depth-velocity  percent  fre¬ 
quencies  are  presented  for  four  channel  categories  (divided,  narrow,  transi¬ 
tional,  and  wide)  in  a  median  year.  Coefficients  are  unadjusted,  adjusted,  and 
adjusted  and  weighted  by  reach. 
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B1 


Tabfe  B1 

Gavins  Point  Canberra  Coefficients  of  Historicai/Operationai  Depth-Veiocity 
Percent  Frequencies  Unadjusted  for  Divided  Channei  Category  in  a  Median 
Year 


OBS  JAN  FEB  MAR  APR  MAY  JUN  JUL  AUG  SEP  OCT  NOV  DEC 


1  0.22230  0.18001  0.20648  0.19410  0.21508  0.21508 

2  0.25125  0.24568  0.21061  0.21661  0.24203  0.24203 

3  0.28059  0.28035  0.25574  0.24138  0.25555  0.25555 

4  0.25387  0.27904  0.27672  0.28194  0.29264  0.29264 

5  0.26133  0.29340  0.32908  0.32810  0.33409  0.33409 

6  0.33307  0.37245  0.39204  0.41350  0.44133  0.44133 

7  0.42286  0.44586  0.48074  0.50149  0.56624  0.56624 

8  0.37470  0.40292  0.42981  0.45008  0.52829  0.52829 

9  0.38896  0.38523  0.49846  0.52675  0.54131  0.54131 

10  0.39551  0.40710  0.48583  0.52132  0.50313  0.50313 

11  0.36102  0.36249  0.48961  0.53649  0.53348  0.53348 

12  0.42412  0.40670  0.53535  0.53897  0.53714  0.53714 

13  0.39476  0.39426  0.46223  0.45297  0.46203  0.46203 

14  0.35145  0.33037  0.41991  0.42080  0.43828  0.43828 


0.21428  0.17972  0.17770  0.17837  0.18001  0.22309 
0.22712  0.23147  0.23917  0.23609  0.24568  0.24530 
0.26361  0.24980  0.26980  0.26292  0.28035  0.26148 
0.29541  0.27976  0.27789  0.27916  0.27904  0.25114 
0.33975  0.29026  0.29056  0.29049  0.29340  0.26478 
0.44021  0.39297  0.37807  0.38280  0.37245  0.32916 
0.53681  0.47065  0.44888  0.45513  0.44586  0.41912 
0.47448  0.44396  0.41151  0.42195  0.40292  0.36276 
0.53870  0.44950  0.40182  0.41742  0.38523  0.38421 
0.51770  0.45570  0.42197  0.43254  0.40710  0.39298 
0.53070  0.39267  0.37011  0.37765  0.36249  0.35937 
0.53858  0.45672  0.41856  0.42953  0.40670  0.42876 
0.46182  0.41553  0.40212  0.41026  0.39426  0.39003 
0.42145  0.34237  0.33347  0.33632  0.33037  0.34636 


OBS  JAN  FEB  MAR  APR  MAY  JUN 

1  0.14001  0.10891  0.07684  0.07670  0.07202  0.07202 

2  0.15825  0.14864  0.07838  0.08559  0.08105  0.08105 

3  0.19873  0.19073  0.10702  0.10726  0.09623  0.09623 

4  0.21690  0.22900  0.13969  0.15112  0.13294  0.13294 

5  0.22327  0.24079  0.16612  0.17587  0.15177  0.15177 

6  0.29983  0.32206  0.20852  0.23353  0.21123  0.21123 

7  0.40305  0.44338  0.29735  0.32936  0.31517  0.31517 

8  0.33006  0.37270  0.28421  0.31601  0.31435  0.31435 

9  0.32361  0.33824  0.34400  0.38600  0.33617  0.33617 

10  0.28193  0.31085  0.36949  0.42098  0.34433  0.34433 

11  0.23713  0.25729  0.38856  0.45208  0.38098  0.38998 

12  0.25309  0.26520  0.44387  0.47449  0.40076  0.40076 

13  0.23557  0.25709  0.38324  0.39878  0.34472  0.34472 

14  0.20972  0.21542  0.34816  0.37046  0.32700  0.32700 


JUL  AUG  SEP  OCT  NOV  DEC 

0.07697  0.10448  0.10664  0.10619  0.10891  0.14746 
0.08158  0.13457  0.14353  0.14055  0.14864  0.16213 
0.10647  0.16331  0.18207  0.17601  0.19073  0.19434 
0.14394  0.22062  0.22622  0.22544  0.22900  0.22517 
0.16554  0.22891  0.23654  0.23459  0.24079  0.23740 
0.22599  0.32653  0.32428  0.32571  0.32206  0.31096 
0.32047  0.45478  0.44773  0.45034  0.44338  0.37781 
0.30283  0.42930  0.38421  0.39755  0.37270  0.29948 
0.35883  0.41438  0.35645  0.37400  0.33824  0.29747 
0.38002  0.36998  0.32642  0.33883  0.31085  0.25513 
0.40651  0.29850  0.26656  0.27584  0.25729  0.21219 
0.43100  0.32187  0.27748  0.28935  0.26520  0.22613 
0.36957  0.29284  0.26659  0.27636  0.25709  0.20570 
0.33727  0.24128  0.22108  0.22656  0.21542  0.18267 


OBS  JAN 


FEB  MAR 


APR  MAY  JUN  JUL  AUG  SEP  OCT  NOV  DEC 


1 

0.019602 

0.015247 

2 

0.022155 

0.020809 

3 

0.027822 

0.026702 

4 

0.030367 

0.032060 

5 

0.031258 

0.033710 

6 

0.041976 

0.045089 

7 

0.056427 

0.062073 

8 

0.046208 

0.052178 

9 

0.045305 

0.047354 

10 

0.039471 

0.043519 

11 

0.033198 

0.036020 

12 

0.035433 

0.037128 

13 

0.032980 

0.035992 

14 

0.029361 

0.030159 

0.010758  0.010738  0.010083 
0.010973  0.011983  0.011347 
0.014983  0.015016  0.013472 
0.019557  0.021157  0.018611 
0.023257  0.024621  0.021247 
0.029193  0.032694  0.029573 
0.041629  0.046110  0.044123 
0.039790  0.044242  0.044009 
0.048161  0.054040  0.047063 
0.051729  0.058938  0.048206 
0.054399  0.063292  0.053338 
0.062141  0.066429  0.056106 
0.053654  0.055829  0.048260 
0.048742  0.051865  0.045780 


0.010083  0.010775  0.014628 
0.011347  0.011421  0.018840 
0.013472  0.014906  0.022863 
0.018611  0.020151  0.030887 
0.021247  0.023176  0.032047 
0,029573  0.031639  0.045714 
0.044123  0.044866  0.063669 
0.044009  0.042396  0.060102 
0.047063  0.050237  0.058013 
0.048206  0.053203  0.051797 
0.053338  0.056911  0.041791 
0.056106  0.060340  0.045061 
0.048260  0.051740  0.040998 
0.045780  0.047217  0.033780 


2 


0.014930  0.014867  0.015247  0.020644 
0.020094  0.019677  0.020809  0.022699 
0.025490  0.024641  0.026702  0.027208 
0.031671  0.031562  0.032060  0.031524 
0.033115  0.032842  0.033710  0.033236 
0.045399  0.045600  0.045089  0.043534 
0.062683  0.063047  0.062073  0.052893 
0.053789  0.055657  0.052178  0.041927 
0.049903  0.052360  0.047354  0.041646 
0.045699  0.047436  0.043519  0.035718 
0.037318  0.038618  0.036020  0.029707 
0.038848  0.040509  0.037128  0.031658 
0.037322  0.038691  0.035992  0.028798 
0.030951  0.031718  0.030159  0.025574 
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Table  B1  (Continued) 


OBS  JAN  FEB  MAR  APR  MAY  JUN  JUL  AUG  SEP  OCT  NOV  DEC 


1  0.28083 

2  0.31725 

3  0.26979 

4  0.28956 

5  0.33743 

6  0.30070 

7  0.23785 

8  0.20900 

9  0.24670 

10  0.30093 

11  0.25419 

12  0.21342 

13  0.22850 

14  0.19302 


0.31653 

0.34124 

0.32990 

0.36119 

0.37199 

0.33868 

0.23311 

0.21294 

0.21968 

0.21171 

0.22459 

0.22423 

0.21980 

0.18765 


0.21726 

0.25468 

0.27223 

0.30760 

0.30846 

0.28522 

0.21433 

0.21942 

0.23240 

0.25501 

0.24681 

0.23417 

0.22883 

0.20697 


0.22973 

0.26906 

0.28206 

0.30281 

0.32380 

0.28533 

0.21847 

0.23822 

0.24114 

0.25911 

0.25477 

0.24118 

0.23857 

0.22689 


0.22973 

0.26906 

0.28206 

0.30281 

0.32380 

0.28533 

0.21847 

0.23822 

0.24114 

0.25911 

0.25477 

0.24118 

0.23857 

0.22689 


0.19026 

0.22331 

0.23652 

0.26968 

0.27871 

0.25286 

0.24505 

0.21510 

0.20271 

0.23851 

0.23189 

0.20023 

0.21094 

0.18776 


0.22889 

0.26278 

0.28167 

0.31662 

0.31195 

0.29468 

0.23990 

0.23184 

0.23961 

0.26012 

0.26500 

0.23039 

0.23741 

0.20543 


0.32386 

0.34600 

0.33450 

0.36502 

0.38609 

0.34368 

0.22459 

0.21443 

0.22748 

0.21416 

0.22802 

0.21739 

0.22897 

0.21347 


0.31148 

0.33645 

0.32387 

0.35217 

0.36591 

0.33161 

0.22937 

0.20978 

0.21817 

0.21002 

0.22281 

0.22043 

0.22164 

0.18933 


0.31466 

0.33877 

0.32295 

0.35140 

0.36842 

0.33182 

0.22704 

0.20992 

0.21807 

0.21056 

0.22266 

0.21798 

0.22364 

0.19329 


0.31653 

0.34124 

0.32990 

0.36119 

0.37199 

0.33868 

0.23311 

0.21294 

0.21968 

0.21171 

0.22459 

0.22423 

0.21980 

0.18765 


0.29985 

0.34192 

0.30922 

0.32253 

0.35469 

0.32312 

0.22950 

0.21071 

0.24815 

0.26174 

0.25279 

0.20614 

0.21687 

0.18375 


OBS  JAN  FEB  MAR  APR  MAY  JUN  JUL  AUG  SEP  OCT  NOV  DEC 


1  0.13953  0. 

2  0.17020  0. 

3  0.15559  0. 

4  0.18180  0. 

5  0.21988  0. 

6  0.22645  0. 

7  0.18536  0. 

8  0.16653  0. 

9  0.20407  0. 

10  0.25970  0. 

11  0.23236  0. 

12  0.21266  0. 

13  0.21526  0. 

14  0.17418  0. 
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.14944  0.( 
.17812  0.1 
.19040  0.1 
20034  0.1 
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.16126  0.1 
17550  0.1 
18265  0.1 
16178  0.1 


\  0.07025 
I  0.08884 
:  0.10011 
i  0.11701 
0.12986 
I  0.13225 
'  0.10479 
0.11683 
0.12276 
i  0.13763 
i  0.14334 
0.14791 
10.15534 
0.15327 


0.05788  0.( 
0.07336  O.C 
0.08352  O.C 
0.10367  0.1 
0.11120  0.1 
0.11659  0.1 
0.11693  0.1 
0.10494  0.1 
0.10267  0.1 
0.12603  0.1 
0.12979  0.1 
0.12216  0.1 
0.13664  0.1 
0.12618  0.1 


I  0.06999  0. 
)  0.08676  0. 
’  0.09997  0. 
r  0.12234  0. 
10.12510  0. 
I  0.13657  0. 
I  0.11506  0. 
i  0.11369  0. 

'  0.12198  0. 
1  0.13815  0. 
)  0.14909  0. 
)  0.14128  0. 
i  0.15458  0. 
\  0.13877  0. 


.11888  0. 
.13714  0. 
.14251  0. 
.16931  0. 
.18586  0. 
.19120  0. 
,12930  0. 
12623  0, 
13901  0. 
.13654  0. 
.15399  0. 
.16003  0. 
.17896  0. 
.17310  0. 


.12004  0. 
.14001  0. 
.14488  0. 
.17151  0. 
.18495  0. 
.19371  0. 
,13865  0. 
12966  0. 
13998  0. 
.14059  0. 
.15799  0. 
.17038  0. 
.18189  0. 
,16120  0. 


.11977  0, 
13924  0. 
14268  0. 
16902  0. 
18392  0. 
19143  0. 
13555  0. 
12814  0. 
13819  0. 
13921  0. 
15593  0. 
16639  0. 
18126  0. 
16253  0. 


M2353  0.15148 
.14380  0.18651 
.14944  0.18131 
.17812  0.20588 
.19040  0.23499 
.20034  0.24741 
.14270  0.18185 
.13327  0.17071 
.14273  0.20870 
.14351  0.22967 
.16126  0.23495 
.17550  0.20344 
.18265  0.20045 
.16178  0.16243 


OBS  JAN 

1  0.030697 

2  0.037445 

3  0.034229 

4  0.039995 

5  0.048373 

6  0.049820 

7  0.040780 

8  0.036637 

9  0.044895 

10  0.057134 

11  0.051120 

12  0.046785 

13  0.047356 

14  0.038319 


0.027176 

0.031636 

0.032877 

0.039186 

0.041887 

0.044074 

0.031393 

0.029320 

0.031401 

0.031572 

0.035478 

0.038609 

0.040183 

0.035592 


0.014616 

0.018501 

0.021257 

0.026149 

0.027216 

0.029082 

0.022616 

0.023673 

0.026028 

0.029797 

0.030549 

0.031593 

0.032780 

0.030759 


0.015455 

0.019545 

0.022025 

0.025741 

0.028569 

0.029095 

0.023053 

0.025702 

0.027008 

0.030278 

0.031534 

0.032540 

0.034175 

0.033720 


0.012734 

0.016138 

0.018374 

0.022807 

0.024464 

0.025651 

0.025725 

0.023087 

0.022586 

0.027727 

0.028554 

0.026876 

0.030061 

0.027760 


0.012734 

0.016138 

0.018374 

0.022807 

0.024464 

0.025651 

0.025725 

0.023087 

0.022586 

0.027727 

0.028554 

0.026876 

0.030061 

0.027760 


0.015398 

0.019088 

0.021993 

0.026914 

0.027522 

0.030046 

0.025313 

0.025012 

0.026835 

0.030394 

0.032799 

0.031082 

0.034007 

0.030529 


0.026153 

0.030170 

0.031353 

0.037247 

0.040890 

0.042065 

0.028447 

0.027770 

0.030582 

0.030038 

0.033878 

0.035206 

0.039372 

0.038082 


0.026410 

0.030803 

0.031873 

0.037732 

0.040689 

0.042616 

0.030504 

0.028525 

0.030796 

0.030930 

0.034757 

0.037483 

0.040015 

0.035464 


0.026349 

0.030632 

0.031389 

0.037184 

0.040461 

0.042115 

0.029821 

0.028191 

0.030401 

0.030627 

0.034305 

0.036607 

0.039877 

0.035757 


0.027176 

0.031636 

0.032877 

0.039186 

0.041887 

0.044074 

0.031393 

0.029320 

0.031401 

0.031572 

0.035478 

0.038609 

0.040183 

0.035592 


0.033325 

0.041033 

0.039889 

0.045295 

0.051699 

0.054430 

0.040007 

0.037556 

0.045914 

0.050526 

0.051689 

0.044756 
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Appendix  B  Gavins  Point  Canberra  Coefficients 


Table  B1  (Concluded) 


OBS  JAN  FEB  MAR  APR  MAY  JUN  JUL  AUG  SEP  OCT  NOV  DEC 

1  0.29929  0.31224  0.36315  0.36918  0.35833  0.35833  0.39465  0.30581  0.30941  0.30794  0.31224  0.29031 

2  0.31180  0.40083  0.43304  0.44126  0.40143  0.40143  0.46186  0.37828  0.39258  0.38834  0.40083  0.30587 

3  0.30568  0.37020  0.48268  0.48770  0.48129  0.48129  0.51708  0.36007  0.37520  0.37608  0.37020  0.28165 

4  0.33732  0.42884  0.54054  0.52578  0.47669  0.47699  0.55865  0.41393  0.43496  0.43311  0.42884  0.31737 

5  0.31716  0.39683  0.50729  0.50925  0.52474  0.52474  0.51517  0.40613  0.40589  0.41044  0.39683  0.29028 

6  0.33499  0.42813  0.56063  0.57068  0.56094  0.56094  0.54044  0.44840  0.44071  0.44593  0.42813  0.31702 

7  0.25627  0.32509  0.49964  0.50413  0.50488  0.50488  0.49435  0.33682  0.33329  0.33897  0.32509  0.25156 

8  0.25612  0.30868  0.48575  0.49020  0.50056  0.50056  0.47612  0.32578  0.31712  0.32389  0.30868  0.24148 

9  0.26723  0.31739  0.48431  0.49218  0.46559  0.46559  0.47989  0.33612  0.32421  0.32803  0.31739  0.25878 

10  0.26097  0.30627  0.47325  0.48872  0.44519  0.44519  0.46274  0.32873  0.31183  0.31688  0.30627  0.26553 

11  0.26740  0.31854  0.51388  0.53064  0.46043  0.46043  0.49578  0.34026  0.32850  0.33351  0.31854  0.26715 

12  0.28862  0.30875  0.53042  0.52948  0.45574  0.45574  0.51920  0.32469  0.31471  0.31974  0.30875  0.29632 

13  0.29307  0.31111  0.51718  0.50481  0.52668  0.52148  0.52148  0.33130  0.31877  0.32508  0.31111  0.29442 

14  0.30515  0.31011  0.47760  0.47630  0.50447  0.50447  0.49167  0.32528  0.31295  0.31507  0.31011  0.29757 


OBS  JAN  FEB  MAR  APR  MAY  JUN  JUL  AUG  SEP  OCT  NOV  DEC 

1  0.22864  0.19575  0.13058  0.13377  0.12295  0.12295  0.13993  0.17247  0.18974  0.18480  0.19575  0.22293 

2  0.24256  0.25590  0.15856  0.16282  0.14026  0.14026  0.16676  0.21724  0.24515  0.23732  0.25590  0.23918 

3  0.25081  0.24927  0.18641  0.18980  0.17736  0.17736  0.19691  0.21809  0.24712  0.24240  0.24927  0.23229 

4  0.28003  0.29216  0.21121  0.20703  0.17773  0.17773  0.21524  0.25367  0.28985  0.28245  0.29216  0.26483 

5  0.27585  0.28325  0.20768  0.21009  0.20498  0.20498  0.20796  0.26076  0.28338  0.28043  0.28325  0.25378 

6  0.30027  0.31493  0.23654  0.24263  0.22582  0.22582  0.22483  0.29671  0.31710  0.31400  0.31493  0.28563 

7  0.20941  0.31558  0.27818  0.28285  0.26822  0.26822  0.27139  0.29411  0.31646  0.31498  0.31558  0.20402 

8  0.14933  0.25842  0.32397  0.32946  0.31855  0.31855  0.31311  0.31079  0.27354  0.28726  0.25842  0.13903 

9  0.09873  0.21008  0.37170  0.38066  0.34097  0.34097  0.36316  0.26765  0.22407  0.23588  0.21008  0.09343 

10  0.06598  0.17339  0.38920  0.40502  0.34935  0.34935  0.37523  0.23346  0.18631  0.19883  0.17339  0.06473 

11  0.01128  0.10323  0.49396  0.51400  0.42231  0.42231  0.46990  0.16756  0.11849  0.13198  0.10323  0.01409 

12  0.02821  0.08598  0.52373  0.52683  0.42938  0.42938  0.50549  0.14657  0.09948  0.11258  0.08598  0.03218 

13  0.05550  0.06324  0.50140  0.48539  0.51789  0.51789  0.51307  0.12715  0.07731  0.09106  0.06324  0.05908 

14  0.07083  0.05216  0.45341  0.44832  0.50319  0.50319  0.47399  0.11457  0.06516  0.07767  0.05216  0.07251 


OBS 

JAN 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

OCT 

NOV 

DEC 

0.075451 

0.06460 

0.04309 

0.04414 

0.04057 

0.04057 

0.04618 

0.05691 

0.06261 

0.06098 

0.06460 

0.073567 

2 

0.080046 

0.08445 

0.05233 

0.05373 

0.04629 

0.04629 

0.05503 

0.07169 

0.08090 

0.07832 

0.08445 

0.078929 

3 

0.082767 

0.08226 

0.06152 

0.06263 

0.05853 

0.05853 

0.06498 

0.07197 

0.08155 

0.07999 

0.08226 

0.076656 

4 

0.092408 

0.09641 

0.06970 

0.06832 

0.05865 

0.05865 

0.07103 

0.08371 

0.09565 

0.09321 

0.09641 

0.087393 

5 

0.091030 

0.09347 

0.06853 

0.06933 

0.06764 

0.06764 

0.06863 

0.08605 

0.09352 

0.09254 

0.09347 

0.083746 

6 

0.099089 

0.10393 

0.07806 

0.08007 

0.07452 

0.07452 

0.07419 

0.09791 

0.10464 

0.10362 

0.10393 

0.094257 

7 

0.069106 

0.10414 

0.09180 

0.09334 

0.08851 

0.08851 

0.08956 

0.09706 

0.10443 

0.10394 

0.10414 

0.067327 

8 

0.049280 

0.08528 

0.10691 

0.10872 

0.10512 

0.10512 

0.10333 

0.10256 

0.09027 

0.09479 

0.08528 

0.045879 

9 

0.032581 

0.06933 

0.12266 

0.12562 

0.11252 

0.11252 

0.11984 

0.08833 

0.07394 

0.07784 

0.06933 

0.030830 

10 

0.021772 

0.05722 

0.12844 

0.13366 

0.11529 

0.11529 

0.12383 

0.07704 

0.06148 

0.06561 

0.05722 

0.021361 

11 

0.003721 

0.03407 

0.16301 

0.16962 

0.13936 

0.13936 

0.15507 

0.05529 

0.03910 

0.04355 

0.03407 

0.004649 

12 

0.009309 

0.02837 

0.17283 

0.17386 

0.14170 

0.14170 

0.16681 

0.04837 

0.03283 

0.03715 

0.02837 

0.010618 

13 

0.018314 

0.02087 

0.16546 

0.16018 

0.17090 

0.17090 

0.16931 

0.04196 

0.02551 

0.03005 

0.02087 

0.019498 

14 

0.023374 

0.01721 

0.14963 

0.14795 

0.16605 

0.16605 

0.15642 

0.03781 

0.02150 

0.02563 

0.01721 

0.023928 
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